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.SUMMARY 


Tht*  t •<  tin*  Si  intiiimn.il  'l  et  In. h  ul  Summary  of  tin*  Si-mmn  UiMi-imhiunott  «  roup 

•  if  l.ltn  idii  l.iil  iir.itiTV.  It  ill* fit* i'll ii*.i  ri^M'tirrlt  iii»mi*  tiiti*il  wtllt  tin*  tliinrv  :nnl  |ii‘iii  1 1 1'**  i*f  il> * 1 1*1* 1 1 •  >n. 
Ini  itlmt  mill  tili-m i f n-fit ton  iif  i*artlii|iiiiki**i  mill  uiul**i  ground  explosion*.  1‘artiridnr  emphasis  is 
|il<n*i*tl  oil  tin*  ruli*  that  digital  tl.it it  processing  au*l  large  army  t i*i- Inn i|in* s  plav  III  such  tmrk, 

A  stildv  lias  I . it  uiinli*  of  some  unusual  i*vi*nts  itt  Asia  which  hill,  hi  ulmosl  fail,  tin*  M  in, 

s  ll 

ilisi  riliiillant,  Tlirsi*  1*1  ills  i  it  tit  i*  a  small  lint  nut  negligible  1 1  ri  >|  >i  >  rt  i  in  iif  events  in  that  region,  Clin* 
t'vwiu  lias  t'iinsi*il  |iartii'iilar  ililfirnltv,  with  ;m  apparent  ili*|>th  |  liasi*  lint  i*vi*ry  utlu*r  as|ii*rt  of  an 
explosion,  A  catalogue  of  |iri*snmi,il  explosions  tit  tin*  Soviet  I  iiion  away  fniir  tin*  regular  tost 
situs  is  |in*si*nti*il.  A  time  iliiniaiu  analysis  of  seismograms  fiimi  a  largo  Siberian  i*;t  rt  In  |ii;i  k<* 
lias  ln*i*it  ininli'.  I'liis  i*vi*nt  lias  special  interest,  being  in  a  region  of  only  inilil  scismirit v.  Cor* 
tain  fuatii res  in  (lit*  deconvolved  record  suggest  that  *npiuro  dimensions  ami  velocity  ran  In*  ox- 
trartrd  from  the  suitr  of  rreords.  Tin*  (|iii*stion  of  tin*  cause  of  discriinimition  is  reopened  with 
a  discussion  of  conclusions  that  can  In*  drawn  from  curth<|iinkrs  of  apparently  very  small  dimen¬ 
sions  which  discriminate  with  east*.  We  conclude  that  the  partition  of  energy  between  I’-wavos 
ami  Rayleigh  waves  may  Ik*  tin*  central  cause,  A  study  has.  been  completed  of  R  phases  from 
cratering  and  fully  contained  explosions  in  Nevada.  We  found  that  explosions  that  break  surface 
are  richer  in  R^  than  fully  contained  events, 

Kxtensive  work  has  boon  done  at  Lincoln  on  the  multipathing  of  Rayleigh  waves,  and  now  a 
ray-tracing  program  has  been  used  on  a  globe  with  oontinont/oooan  distinctions.  The  results 
show  striking  agreement  with  observation,  and  open  up  the  possibility  that  wc*  can  more  adequately 
predict  expected  surface  wave  magnitudes  from  hypothetical  sources  on  a  real  earth.  Travel¬ 
time  studies  up  to  the  present  huve  hud  to  make  assumptions  about  neat*  source  heterogeneities 
which  may  not  be  valid.  Wc  have  therefore  conducted  a  P-wave  travel -time  study  which  use 
only  deep-focus  earthquakes.  A  rather  surprising  result  emerges.  Rays  which  huve  penetrated 
the  deepest  mantle  show  very  systematic  deviations  from  standard  travel-time  tables  of  up  to 
±1  sec.  This  result  is  of  interest  in  characterizing  regions  of  the  deep  mantle  and  also  in  alleinpls 
to  predict  the  value  of  seismic  observatories  at  distances  beyond  90r.  A  cognate  result  has  been 
found  from  recent  array  studies  of  dT/dA  and  azimuth  anomalies.  It  appears  inevitable  that  there 
are  large-  and  small-scale  lateral  variations  of  significant  amplitude  throughout  the  whole  mantle, 
but  particularly  near  the  core- mantle  boundary. 

Work  continues  on  describing  the  structure  beneath  l.ASA,  both  in  terms  of  a  Chernov-type 
approach  in  which  the  amplitude  and  phase  fluctuations  observed  at  individual  seismometers  are 
related  by  a  model  in  which  the  earth  beneath  l.ASA  is  a  randomly  varying  medium  with  small 
variations  in  refractive  index,  and  also  by  an  analysis  of  S-waves  generated  near  the  surface  by 
interfaces  such  as  the  sediment/basemont  boundary.  A  routine  technique  is  proposed  for  quality 
control  of  l.ASA,  or  any  other  digital  seismic  data,  by  use  of  pseudo-random  input  to  tin*  cali¬ 
bration  coil  and  spectral  analysis  of  the  resulting  output  tape. 

The  increasing  availability  of  NORSAR  data  and  the  interchange  of  scientists  with  the  NORSAR 
data  center  lead  to  an  increasing  amount  of  research  into  NORSAR  capabilities.  Studies  of  NORSAR 
signal  amplitudes  have  concentrated  on  the  degree  of  coherence  in  short  "period  dala  across  the 
array,  and  contrasts  are  made  nth  amplitude  scatter  at  l.ASA.  We  find  that  NORSAR  body  wave 


iv"  •  «»y 


ammtuto*  of  u,  even,*  are  generally  l<ivu>r  kM|{|1  „ 

I.igl, -Impedance  vft|,,e„  *  ‘ ,  '  ,  '  *  f,,r 

N'OHSAH  dot. -i  t ton  processor  is  .  i  .  1  *  fsuH‘”,,,"rin  rnio  m 

term*  of  -irol  characierlsties  of  7"  T,"  *n  *'  "re  considered  in 

«•>'"  nro  con, n, sled  and  „  ,s  . *  *'*  ^ “  f,,r  '*ASA  """  *>■«*'«  -•«'« -p-rlod 

rro.MM.,u  y  signal  i„  favor  of  I  en-rgy’  '  “r  °  *°  f"",r  "n*  »»'«AK  high- 

. — . 

•surface  of  the  cnrih  lias  needed  sonw  nr.  '  '  U  °f  1  H‘S,>  fl"rtll,'*'"nx  on  the 

integrals  which  do  no,  converge  .'  The  eoMocation^of  ru'ra7dV'long-rl,M! l" 

Pomeroy  type  at  Fairbanks  \i,  u.  in  ‘  B  pi  riod  .wersniomeler  of 

Of  instruments  recording  in  somewhat  r<''ntiVt‘  mt,r,,s 

intervals  is  being  examined.  y  “  A  Population  of  overlapping  time 

tance of^sln^I.AS/Vs'deteciio^og^i^c"^^*0!  nCtW°r^  CaP°bUity  ^  th.  impor- 

e,mv  “  i,  “ ;;rnjT  ;;,ih  somo  "r"“"y  **»» 

Lincoln  LnWn.or,  Jac,  Z\o  t  Z  T  “*  ^  Month  ,„P  .„lch 

—  tho  rrr;„  ™„“ved  ,h:  co"“,,or"uon  °r ~~ 

short-period  »*»,  ior  ,  ,  -™nUl  perlo„  ^oT.lll'L  **''«“  ^  ™ 

second  phase,  surface  waves  will  be  analyzed  We  ih  h  lt>b"ec  soon  and,  in  the 

contributions  that  can  be  mode  bv  in  r  a  i  °n  °P°  U>  reach  uscful  conclusions  on  the 

be  made  by  individual  stations  to  a  seismic  monitoring  network. 

D.  Davies 
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GLOSSARY 


AI.I’A 

Alaskan  lamp- Period  Array 

CCS 

United  States  Coast  and  Oeodotic  Survey 

CASA 

l.arpe  Aperture  Seismic  Array 

NO  A  A 

National  Oceanoprapbic  and  Atmospheric  Administration 

NOHSAK 

Norwepian  Seismic  Array 

NOS 

National  Ocean  Survey 

l»l)K 

Preliminary  Determination  of  epicenters 

SAAC 

Seismic  Array  Analysis  Center 

SATS 

Semiannual  Technical  Summary 

WWSSN 

World-Wide  Standard  Seismopraphic  Network 
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SEISMIC  DISCRIMINATION 


I.  THE  SEISMIC  SOURCE 

A.  SOME  INTERESTING  CENTRAL  ASIAN  EVENTS  ON  THE  DIAGRAM 

Se  smlc  discrimination  between  Mrthqw.hr*  and  o«plo*inn*  I*  c«mpllrair.|  by  tie  oMrrvittnn 
of  occasional  event*  with  conflicting  discrimination  parameter*.  In  particular,  ihr  mn»f  -He 
„*od  discriminant,  the  m,,  diagram.  classifies  a  small  number  of  earthqiMLr*  a.  «*Mdn*lon« 
(for  example,  possibly  a  few  of  the  earthquake*  occurring  In  the  laguna  SiU*ti  atvl  Nevada- 
Arizona  region*  of  the  United  State*,1  an  l  ce-talnly  two  of  the  event*  utu-lle.l  herd.  Eurthrr 
complication*  arise  when  the  *urface  wave*  from  an  explo*ion  are  masked  by  the  surface  wave* 
of  a  large  earthquake,  or  when  the  bo.ly  wave*  of  an  explosion  are  perturlw.l  to  re»emble  short- 
period  earthquake  radiation  through  the  detonation  of  several  explosion*  over  a  short  time  spin. 
Here,  the  existence  of  some  Tibetan  and  southern  Husslan  events  that  have  anomalously  low 
surface-wave-to-body-wave  magnilu<le  ratios  are  documented,  Further,  one  of  the  southern 
Russian  events  (event  KA  in  Table  1-1)  is  used  to  illustrate  the  problems  that  arise  when  an  inter¬ 
fering  event  is  present  and  when  multiple  detonations  are  considered  possible.  Mar  sits  II  and 
Basham2  have  previously  determined  the  general  behavior  of  M^ntp  for  Central  Asm.  and  their 
relationship  is  used  for  purposes  of  comparison. 


TABLE  M 

EVENTS  USED  IN  THIS  STUDY  (CGS  LOCAT 

IONS) 

Event 

Origin 

Latitude 

Lonitude 

Depth 

(T). 

M* 

Symbol 

Day 

Time  (GMT) 

(°N) 

(“E) 

b 

s 

KA 

5-1-69 

04:00:08.7 

44.0 

77.9 

53 

4.9 

2.4 

KB 

12-9-69 

13:41:09,0 

40. 1 

70.7 

N 

4  9 

3.6 

KC 

7-1-68 

19:14:54.7 

44.0 

79.3 

N 

4.8 

3.6 

KD 

7-20-69 

04:34:14.9 

39.8 

77.8 

N 

4.9 

4.2 

TA 

7-15-68 

05:09:05.9 

30.3 

95.0 

22 

4.8 

2.9 

TB 

7-14-68 

18:12:40.9 

30.3 

94.8 

22 

4.9 

3.3 

TC 

7-13-68 

05:05:54.2 

30.3 

94.6 

N 

5.0 

3.4 

*  Eac!  M 

value  is  bas 

ed  on  the  results  from  at  least  four  individual  stations. 

Tne  location  of  the  stations  and  events  used  in  this  study  are  plotted  on  a  map  of  Kurasia  in 
Fig.  1-1.  These  data  are  a  subset  of  data  being  used  in  a  network  study  of  Central  Asian  seis¬ 
micity  (SATS,  31  December  1971,  DOC  A  D- 737092);  station  data  are  given  in  that  report.  The 
four  events  starting  with  a  K  are  located  in  the  southeastern  corner  of  Kazakh,  USSR.  They 
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Section  I 


locate  within  10  of  each  other.  The  other  three  events  discussed  are  closely  grouped  in  an  area 
of  Tibet  alnuit  H'Okm  northeast  of  the  capital,  Lhasa.  'They  are  prefixed  by  the  letter  'I'.  All 
?  events  have  approximately  the  same  body  wave  magnitude,  4.9  ±  0,1,  and  are  listed  as  shallow 
focus  ear*'  , ’takes  in  t he  NOS  bulletin  (that  is,  they  are  assigned  focal  depths  between  20  and 
70km  i-  I'M  1 :  catalog).  Kpicentral  data  determined  by  NOS  are  given  in  Table  1-1.  The 
Kazakh  c.ent  he  on  the  border  of  a  seism  ie  Central  Asia.  Plowshare-type  explosioi  s  have  been 

di .  ‘  this  area  (see  event  dated  4 0-21  - 6 6  in  Table  1-3).  The  Soviet  test  site  at  Serni- 

palatinsk  is  abotit  b'  ,o  the  north.  The  Tibetan  events  are  part  of  an  isolated  cluster  of  seismic 
activity. 

M  values  were  calculated  using  the  scale  proposed  bv  Marshall  and  Masha 
s 

M  -  logj  ()(A)  +  P(T)  +  IT  (A) 

...  •  .he  maximum  amplitude  of  the  vertical  eotrponent  of  the  Kayleigh  wave  at  the  period 

T  measur*-1  an  epieentral  distance  A.  P(T)  corrects  for  dispersion;  and  M’(A)  corrects  for 
the  average  effects  of  scattering,  geometrical  spreading,  lateral  refraction  and  anelastic  atten¬ 
uation  of  Itavleiidi  waves,  The  M  values  for  the  7  events  are  included  in  Table  1-1  and  range 

S 

from  2. -I  up  to  *1.2.  Thus,  even  though  m.  varies  by  only  0,2,  Mfi  varies  by  i.R.  Mased  on  the 

regional  seismicity  of  Kirghiz  for  events  with  nij  -  4.0  (see  Kef.  3),  the  *1  events  in  that  region 

would  represent  approximately  10  percent  of  the  total  population  that  would  be  observed  over  a 

comparable  period  of  time.  Figure  1-2  shows  how  the  events  plot  on  an  Mfi:m^  diagram.  The  two 

line  drawn  on  the  graph  denote  the  trends  of  earthquake  and  explosion  population  for  Central 

Asia.4-  Data  used  to  calculate  the  trend  lines  were  obtained  from  the  same  stations  and  analyzed 

in  tuc  same  manner  as  the  data  used  here.  None  of  the  earthquakes  used  in  that  study  with  an 

m  •  4  (>  scattered  to  the  left  of  the  event  KM.  Furthermore,  the  inclusion  of  our  events  would 
b 

not  cause  the  trend  lines  to  be  appreciably  altered.  The  events  located  in  Kazakh  (denoted  by 
KB  and  KC),  along  with  the  Tibetan  events  ('I’M  and  TO,  could  not  be  classified  as  either  arti¬ 
ficial  or  natural  on  the  basis  of  M  ;m,  alone.  The  events  KA  and  TA  clearly  fall  into  the  explo¬ 
sion  population.  The  three  Tibetan  events  are  presumably  earthquakes.  They  occur  in  an  un.- 
likelv  geographic  location  for  explosions,  and  have  the  short-period  complexity  and  Lg:  P  energy 
ratios  characteristic  of  earthquakes.  Mor ’over,  many  similar  events  have  occurred  historically 
within  a  degree  or  so  of  their  epicenters.  The  Kazakh  events,  with  the  exception  of  KA,  have 
short-period  records  typical  of  earthquakes, 

A  difference  in  magnitude  of  l.R  cor  responds  to  an  amplitude  ratio  of  approximately  h0,  or 
ISrIM.  Such  large  amplitude  variations  for  these  events,  having  the  same  body  wave  magnitude, 
are  illustrated  in  Fig.  1-3.  The  records  presented  in  the  figure  arc  at  gains  and  epieentral  dis¬ 
tances  that  make  them  directly  comparable.  As  semi  in  Fig.  1-3,  the  seismograms  are  typical  of 
first-  and  second-zone  continental  Kayleigh  waves  from  shallow  sources.  The  Kavleigh  wave  for 
event  KA  is  immersed  in  a  Kayleigh  wave  from  an  earthquake  in  Tonga.  Figure  1-4  shows  the 
three  LI*  components  for  this  event  as  recorded  at  the  WWSSN  observatory  at  Kabul.  The  pres¬ 
ence  of  the  Kavleigh  wave  from  the  southern  Kussia  event  is  clearly  visible  because  of  its  high- 
f requeue v  nature  relative  to  the  low- freoueney  nature  of  the  interfering  event.  It  is  interesting 
to  note  tliat  the  Kazakh  event  would  not  have  been  extractable  if  the  record  1  in  cl  been  produced  by 
a  LASA  tvpe  long-period  or  a  l.’inmnl  Ogdensburg-type  instrument,  since  (.-sec  energy  (the 
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dominant  period  or  tlie  Kazakh  event)  is  down  in  amplitude  by  more  than  an  order  of  magnitude 
from  40-see  energy  (the  dominant  period  of  the  Tonga  event).  Fortunately,  WWSSN  long-period 
instruments  are  rather  broadband:  6-  and  40-sec  energy  have  the  same  magnification  which  is 
only  a  factor  of  two  less  than  the  peak  magnification  that  occurs  near  20  sec  (see  Hef.  4). 

The  ?  events  listed  as  earthquakes  by  NOS  have  been  shown  to  span  the  earthquake  and  ex¬ 
plosion  populations  on  the  M  :mj  diagram.  The  Tibetan  events  which  we  take  to  be  earthquakes 
lie  nearer  the  explosion  population.  The  Kazakh  event  that  lies  within  100  km  of  the  city  of  Alma 
Ata  anti  which  occurred  on  1  Mav  1969  at  04  00  GMT  would  be  presumed  to  be  an  explosion  on  the 
basis  of  M  :  m^.  However,  the  event  is  travel-time  constrained  to  a  depth  of  53  km  in  the  NOS 
bulletin. 

To  examine  the  behavior  of  trr  time  residuals  as  a  function  of  depth,  thirteen  widespread 
WWSSN  stations  where  first  arrival  times  for  the  event  KA  could  be  precisely  determined  were 
used  to  recompute  a  hypoccnter.  Table  1-2  gives  various  hypoeentral  characteristics  based  on 
these  picks,  liven  though  the  minimum  standard  deviation  (as  determined  by  our  location  nro- 
gram)  occurs  for  a  depth  of  I'M  km  at  04:00:06  GMT,  the  error  for  a  surface  focus  at  04:00:00 
GMT  is  not  significantly  different.  The  epicenters  corresponding  to  the  various  hvpoeentcrs 
differ  bv  less  than  3  km.  Tims,  it  is  as  likely  that  the  event  occurred  at  the  surface  on  the  hour 
as  at  4  0- km  depth  at  6  see  after  the  hour, 
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TABLE  1-2 

CHARACTERISTICS  FOR  EVENT  KA 

Latitude 

Longitude 

GMT 

Standard 

(°  N) 

(°E) 

Depth 

Time 

Deviation 

44.  10 

77.  84 

0 

04:00:00 

0.61 

44.09 

77.83 

20 

04:00:03 

0.58 

44.09 

77.  82 

39.  1 

04:00:06 

0.57 

44.09 

77.82 

40 

04:00:06 

0.57 

44.08 

77.80 

60 

04:00:08 

0.59 

This  event  exhibits  other  short-period  characteristics  that,  are  more  typical  of  explosions 
than  earthquakes  hut  arc  not  strong  evidence  either  wav.  It  lias  congressional  first  arrivals 
where  first  motion  can  be  clearly  seen.  The  first  cycle  of  motion  has  a  dominant  period  of 
0.8  cc,  and  the  umplitu''  ratio  of  the  first  dilatational  peak  to  the  first  compression;!  1  peak  is 
about  5  to  1.  The  seismograms  also  have  explosion- like  complexity.  One  first-zone  record 
showing  I.g  was  available  at  Kabul.  It  had  an  amplitude  slightly  less  than  the  initial  arrival.  Such 
a  ratio  is  not  definitive  as  to  either  a  natural  or  an  artificial  source. 

Thus  far,  the  sum  of  the  evidence  on  event  KA  points  to  its  being  an  explosion.  However, 
this  event’s  short-period  records  contain  an  apparent  depth  phase.  ‘Figure  1-5  presents  the  short- 
period  records  from  I  ASA  and  the  arrays  Yellowknife  (YKA)  and  Warramunga  (WHA).  The  ar¬ 
rival  at  7. 5  sec  after  1*  appears  to  lie  pP;  it  seems  to  he  inverted  and  to  have  the  same  wave 
shape.  The  arrival  at  I  ASA  at  about  1  1  sec  after  I*  may  be  sP.  pl’-P  time  puts  the  focus  at 
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2n  to  25  km;  sP-pP  time  puts  the  focus  at  35  to  45  km.  With  the  exception  of  this  phase,  event 
KA  should  be  classified  as  an  explosion.  !f  the  event  were  an  explosion,  a  possible  explanation 
of  the  second  arrival  would  be  a  second  event.  Multiple-explosion  experiments  have  been  con¬ 
ducted  in  the  Russian  Plowshare  program.5  It  seems  possible  that  the  second  arrival  could  have 
been  produced  in  such  a  manner.  If  the  event  is  an  earthquake,  it  most  likely  belongs  in  the 
high  stress  drop  -  small  source  dimension  class  of  events  which  occasionally  occur  in  western 
North  \merica  and  possibly  in  Tibet.  In  either  case,  the  problem  of  discrimination  for  Cen¬ 
tral  Asia  is  made  more  difficult  by  the  existence  of  such  onts. 

I  would  like  to  express  gratitude  to  Peter  Marshall  for  a  prepublication  copy  of  the  paper2 
he  co-authored. 

T.  Landers 

B.  SOME  PRESUMED  EXPLOSIONS  IN  THE  SOVIET  UNION 

1'able  1-3  lists  events  which  may  prove  helpful  in  discrimination  studies.  All  events  on  this 
list  are  believed  to  be  genuine,  although  down  at  the  very  low  magnitude  end  only  LASA  has  con¬ 
tributed.  Care  has  been  taken  to  eliminate  PP  phases.  Several  of  the  events  are  indicated  by 
the  M^: m^  discriminant  to  be  explosions,  and  certain  events  can  be  associated  with  experiments 
described  in  the  literature.5'^  Others  are  simply  events  in  aseismic  regions.  The  list  should 
not  be  regarded  either  as  exclusively  containing  explosions  nor  as  a  comprehensive  listing  of 
explosions  away  from  the  two  regular  test  sites  -  Semipalatinsk  and  Novaya  Zemlya. 

I).  Davies 

C.  DECONVOLUTION  AND  THE  SOURCE  OF  A  SIBERIAN  EARTHQUAKE 

The  interpretation  of  seismograms  in  order  to  infer  the  nature  of  the  seismic  source  is  made 
more  confidently  if  the  effects  of  wave  propagation  and  the  recording  instrument  are  considered. 
Current  trends  in  source  studies  apply  corrections  for  these  effects  in  the  frequency  domain,  and 
the  subsequent  interpretation  is  based  on  the  resulting  spectrum.  Although  in  some  respects  the 
spectral  approach  is  preferred,  the  details  of  the  time  history  of  the  source  are  lost  upon  trans¬ 
formation,  e.g.,  the  earthquahe  could  either  begin  or  end  with  high-frequency  motion  and  yield 
the  same  amplitude  spectrum  in  either  case.  Propagation  and  recording  effects  can  be  removed 
by  well-established  time-domain  techniques,  and  in  what  follows  we  describe  an  attempt  to  apply 
these  techniques  to  a  specific  Asian  earthquake. 

If  P(t)  represents  the  ground  motion  at  a  station  due  to  an  incident  compressional  (P)  wave, 
the  seismogram  S(t)  can  be  represented  symbolically  as  a  convolution 

Pit)  *  lit)  =  S(t)  (F  j, 

where  I(t)  is  the  displacement  impulse  response  of  the  recording  instrument.  The  method  of  de- 
convolution,  or  inverting  Eq.  (1-1)  in  order  to  gain  P(t),  is  well  ..nown.7  An  alternate  method 

may  be  used  to  save  time  if  many  seismograms  recorded  on  similar  instruments  are  to  be  de- 
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convolved.  If  the  inverse  of  lit),  1  It),  is  computed  such  that 

I  1  (t)  *  lit)  =  <5(t)  (1-2) 

where  <5 ( t )  is  the  delta  function,  then  P(t)  may  be  computed  by  the  convolution 

P(t)  =  I-1(t)  *  S(t)  .  (1-3) 
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The  advantage  of  the  second  method  is  obvious  -  only  one  inversion  is  required,  and  S(t)  from 
like  instruments  may  be  deconvolved  through  a  convolution  with  1  (t). 

These  techniques  are  demonstrated  in  Figs.  1-6  and  1-7.  Figure  I-6(a)  shows  the  impulse 
response  of  a  long-period  seismometer  used  in  the  WWSSN  -  this  would  be  lit).  .  igure  I-6(b) 
shows  the  function  I_1(t)  computed  by  inverting  Kq.  (1-2),  and  Fig.  1-6(0  shows  the  result  of  a 
convolution  of  the  functions  I(t>  and  l" 1  ( t ).  If  the  inversion  process  were  exact,  Fig.  1-6(0  would 
be  a  delta  function,  and  it  is  seen  that  this  is  not  strictly  the  case.  Although  Fig.  I-6(c)  is  sharply 
peaked,  the  net  area  under  this  curve  s  probably  closer  to  0  than  1  because  Fig.  I-6(b)  is  essen¬ 
tially  a  differentiator;  it  enhances  high  frequency  motion  lost  through  the  instrument  but.  because 
of  its  short  length,  is  not  as  effective  in  correcting  the  long-period  motion.  Figure  l-7(a)  is  the 
trace  of  a  P-wave  recorded  on  a  WWSSN  long-period  instrument  -  this  would  be  S(t).  Figure  l-7(b) 
shows  the  result  of  computing  P(t)  using  Kq.  (1-  1),  i.e.,  a  convolution  of  Fig.  I-6(b)  with  l-7(a). 

In  this  case,  the  filter  1~ 1  (t )  was  25  points  long  and  S(t)  consisted  of  135  samples.  Figure  l-7(e) 
shows  the  function  P(t>  gained  by  inverting  Kq.  (1-1  )  directly.  It  is  clear  that  although  Figs.  I- 7(b) 
and  (c)  are  similar  in  detail,  the  latter  has  retained  relatively  more  long-period  motion. 

At  present,  we  are  using  these  techniques  to  deconvolve  P-waves  recorded  on  long-period 
instruments  of  the  WWSSN  in  an  attempt  to  gain  some  insight  into  active  faulting  processes.  Ue 
magnify  the  70-mm  films  of  the  original  WWSSN  records  by  successively  using  a  viewer-printer 
and  a  conying  machine  charged  with  transparent  sheets.  The  transparent  sheets  are  then  in¬ 
serted  into  the  viewer-printer  and.  finally,  the  phase  of  interest  is  digitized.  The  entire  process 
before  digitization  yields  a  gain  of  about  30  from  the  original  paper  records.  A  digitization  rate 
of  6  samples/sec  is  realizable;  however,  half  that  rate  is  found  adequate  for  routine  use.  Suit¬ 
able  instrument  constants  are  gained  by  matching  the  Heaviside  acceleration  response  given  by 

the  calibration  pulses  found  on  the  original  records. 

Figure  1-8  shows  the  location  and  fault  plane  solution  of  a  large  earthquake  that  occurred  in 
eastern  Siberia  on  18  May  1971.  The  Karthquakc  I  fata  Kcport  gives  the  following  parameters  for 
this  event: 

Origin  time  22:44:4  3.8 

latitude  6  3.95rN 

Ijongitude  146. 11*  K 

Depth  Normal 

Magnitudes  (stations)  m()  =  5.8(26),  Mg  6.6(9) 

In  Fig.  1-8,  the  base  map  and  "main  deep  faults"  are  due  to  Yanshin.8  The  fault  plane  solution, 
based  on  long-period  first-motion  polarities,  is  consistent  with  a  left-lateral,  vertical  strike 
slip  fault  striking  at  an  azimuth  of  325*.  Although  this  azimuth  is  roughly  parallel  to  the  fault 
traces  given  on  the  map,  the  latter  shows  relative  motion  on  nearby  faults  as  right-lateral.  In 
spite  of  this  discrepancy,  we  take  the  northwest-southeast  plane  to  represent  the  faulting  surface. 

In  Fig.  1-9,  we  show  the  deconvolved  vertical  long-period  V  records,  P(t),  computed  using 
the  second  method  [  Kqs.  (1-2)  and  (1-3)]  described  above.  The  traces  from  28  stations  are  ar¬ 
ranged  at  the  appropriate  azimuth  relative  to  the  fault  plane.  These  traces  represent  ground 
motion  at  the  station  with  no  corrections  being  made  for  path  propagation.  However,  differences 
in  station  gains  have  been  accounted  for  so  that  the  relative  amplitudes  are  correct.  The  reader 
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TABLE  1-3 

EVENTS  IN  ASEISMIC  REGIONS  AND  PRESUMED  EXPLOSIONS  IN  SEISMIC  REGIONS  OF  THE  USSR, 
OTHER  THAN  PRESUMED  EXPLOSIONS  IN  EASTERN  KAZAKH  AND  NOVAYA  ZEMLYA  1966-1971 
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TABLE  1-3  (Continued) 


Section  I 


<N  4) 

.  §  R 

C  -1 

.2  «r  ~ 

S  "  t 

—  4)  _Q 
O-  U  C 
X  UO  5 

c  ^ 
O  t  Q 


O  *°  -5 
«  c  e 
■o-o 

-  £  'O 


/o  5 
<  U  -Q 


*-8 

*«  *- 
0)  ^ 

o  ° 

.y  x 
*o  t: 
c  o 
—  c 

p3  E  C 
F.  -x  4> 

5^5 


£  o  § 

O'  >  a- 


t)  Ml  *■ 

-2-0  0 
*5  ^  V 


8  8  8 
z  z  z 


4>  IS  Q£ 

8  2:  .. 

_  Ml 


-b'S  -C  -£ 
S  >  * 


U  —  U  — 7 

U.  M  £ 


l/l  ^  l/l 

o  So 

u  J  u 


O  < 

s 

-IZ  U 


on  go  lo 

o  o  o 

u  u  u 


a  E  (t\ 

4)  .*  Zr 

O  '  •"  o 


- 

o 

>o 

St 

CO 

00 

WO 

— 

CO 

wo 

wo 

o' 

V 

w> 

wo 

cd 

st 

W) 

WO* 

CO* 

o* 

OG 

o 

OG 

o 

OG 

o 

z 

CO 

>o 

o 

o 

UJ 

0 

00 

UJ 

0 

wo 

UJ 

0 

wo 

UJ 

o 

CN 

UJ 

o 

UJ 

0 

UJ 

0 

wo 

3 

z 

s' 

z 

O* 

wo 

Z 

UJ 

0 

O' 

>o 

WO 

WO 

Z 

UJ 

0 

st 

w> 

^1“ 

z 

IS 

z 

s 

z 

UJ 

0 

9 

UJ 

0 

00 

00 

CO 

St 

St 

CN 

wo 

CO 

5 

Z 

0 

rs 

o 

0 

CN 

3 

z 

0 

CN 

WO 

0 

00 

IS 

wo 

0 

•o 

3 

0 

O 

WO 

z 

0 

St 

Z 

0 

00 

o 

IS 

wo 

IS 

V. 

IS 

'O 

wo 

CO 

wo 

O' 

wo 

CN 

CN 

O'* 

wo 

o 

9. 

O 

9 

CN 

CK 

wo 

O'1 

wo 

O' 

wo 

CN 

st 

O'* 

o 

IS 

o 

IS 

o 

WO 

O 

'O 

o 

SO 

'O 

JO 

R 

R 

R 

R 

O' 

O' 

O' 

O' 

rs 

O' 

V 

.a 

4) 

_Q 

bei 

•“ 

E 

E 

E 

4> 

4) 

> 

4) 

4) 

C 

3 

“I 

O 

Z 

U 

4) 

Q 

U 

4) 

a 

WO 

O' 

CN 

CO 

CN 

CN 

1 — 

CN 

5  8  §  § 


rs 

£  -  K 

IS  O 

i  i  - 

_o  |r. 


s  "  I !  J  I 

2s  >-  «  o  ■*- 


—  Q.O  t; 

-  -  6  0 

r-  O  CN 

CN  Tf  (N 


NORSAR  1  July  1968. 


Section  I 


must  remember  tliat  since  epicentral  distances  here  are  in  the  teleseismic  range  (28°  <  A  <  89°  ), 
wo  are  actually  sampling  the  eompressional  motion  from  a  rather  restricted  solid  angle  of  the 
lower  focal  hemisphere.  In  spite  of  these  re  rictions,  there  are  two  features  of  interest  evident 
in  Fig.  1-9.  The  first  is  that  the  initial  motion  is  usually  emergent,  of  low  amplitude  and  long 
period  relative  to  a  sharp,  high-frequency  motion  arriving  later.  It  is  not  clear  to  us  yet  whether 
this  high-amplitude  motion  represents  a  stopping  phase  of  the  fault  rupture  or  a  breakout  phase 
generated  when  the  faulting  reaches  the  earth's  surface.  The  second  feature  of  interest  is  also 
related  to  this  high-amplitude  motion.  It  may  clearly  be  seen  in  the  right-hand  column  of  traces 
in  Fig.  1-9  that  the  position  of  this  phase  changes  as  a  monotonic  function  of  azimuth.  Although 
this  feature  is  not  as  clear  on  the  traces  on  the  left,  we  take  it  as  an  indication  that  the  source 
of  this  earthquake  propagated  from  the  southeast  to  the  northwest.  In  Fig.  I- 10,  we  compare  the 
deconvolved  trace  at  BRK  with  a  similar  trace  (inverted)  at  ATU.  The  precise  estimation  of  the 
elongation  of  the  trace  at  BKS  with  respect  to  ATU  is  difficult,  but  is  taken  to  be  about  2  to  3  sec. 

If  O  is  tiie  azimuth  to  a  station  measured  clockwise  from  the  direction  of  fault  propagation,  then 
0ATU  *  ^  and  O  =  1 05° .  Using  these  values  of  O,  a  P-phase  velocity  of  16  km/sec,  and 
the  assumption  that  the  late-arriving  high-amplitude  phase  represents  a  stopping  phase  occurring 
at  the  end  of  fault  propagation,  we  would  estimate  the  fault  length  to  be  between  30  to  50  km. 

This  result  is  also  based  on  the  assumption  of  horizontal  propagation  of  rupture  on  a  vertical 
fault.  Estimating  the  velocity  of  rupture  depends  entirely  on  an  accurate  determination  of  the 
duration  of  motion  at  any  given  station.  If,  from  Fig.  1-10,  we  take  the  major  motion  to  endure 
for  8  to  10  sec  at  BRK,  a  rupture  velocity  of  4  to  5  km/sec  is  implied.  This  value  is  higher  than 
the  Rayleigh-wave  velocity,  the  maximum  rupture  velocity  predicted  by  various  crack  propagation 
theories.  I/inger  estimates  of  the  duration  of  motion  will  yield  lower  rupture  velocities. 

J.  Filson 

C.  Frasier 


D.  WHAT  CAUSES  DISCRIMINATION? 

It  has  been  known  for  several  years  that  it  is  possible  to  discriminate  between  explosions 
and  shallo.v  focus  earthquakes  at  least  at  levels  where  signal-to-noise  ratio  causes  no  problems. 
The  Mg:m^  discriminant  has  been  by  far  the  most  reliable  technique  -  in  essence,  explosions 
generate  at  least  an  order-of-magnitude-smaller  Rayleigh  waves  (at  20-st  period)  than  do 
earthquakes  of  the  same  body  wave  magnitude  (at  1-scc  period).  The  reason  why  is  important, 
as  it  may  be  possible  to  predict  the  behavior  of  the  discriminant  to  low  magnitudes.  Yet  the 
cause  of  discrimination  is  poorly  understood.  Several  explanations  have  been  advanced  and  these 
may  be  divided  into  two  categories:  those  that  attempt  to  explain  discrimination  by  spectral  dif¬ 
ferences  between  sources  at  1  and  20  sec,  and  those  that  attempt  to  explain  discrimination  by 
differences  in  partitioning  of  energy  between  body  waves  and  Rayleigh  waves  for  different  sources. 
There  has  been  extensive  discussion  of  source  spectra  -  the  energy  partitioning  problem  has  re¬ 
ceived  less  attention,  probably  because  it  is  rather  difficult  to  describe  an  earthquake  mathemat¬ 
ically.  I  propose  to  present  some  evidence  which,  I  believe,  makes  invalid  any  argument  that 
source  dimensions  or  a  source-time  function  alone  are  the  cause  of  discrimination. 

In  several  previous  SATS,  I  have  remarked  that  core  reflections  such  as  PcP  and  PKiKP 
offer  an  important  view  of  the  seismic  socrce.  There  appear  to  be  several  reasons  for  this.  The 
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core  boundary  seems  to  be  a  clean  reflector  of  seismic  energy.  Seismic  rays  do  not  have  to 
bottom  in  a  vertically  or  late’ ally  heterogeneous  region.  Regions  of  multiple  P-wave  arrivals 
are  avoided.  Plates  and  other  near- source  structure  are  left  rapidly  by  near-vertical  rays. 

S-to-P  conversion  near  source  by  transmission  througn  (on  average)  horizontal  layering  is  slight. 
For  a  good  view  of  core  phases,  the  directional  capability  of  an  array  is  highly  desirable,  hut 
it  is  possible  that  PKP  may  also  prove  valuable  in  some  epicentral  ranges. 

The  striking  thing  to  be  noticed  about  the  use  of  core  phases  is  how  often  they  sharpen  up 
the  short-period  view  of  the  source  when  all  these  effects  mentioned  in  the  last  paragraph  are 
minimized  Figure  I  ll  shows  such  a  case.  The  duration  of  radiation  of  short-period  energy, 
apparently  many  seconds  for  P.  reduces  drastically  for  PcP.  In  fact.  PcP  'rcm  I^ngshot.  with 
a  comparable  body  wave  magnitude,  and  from  a  relatively  nearby  region  lasts  as  long  as.  if  not 
longer  than.  PcP  from  this  earthquake.  It  may  be  argued9  that,  with  certain  provisos,  the  fault 
zone  for  this  earthquake  has  a  maximum  half-dimension  of  1.5  km.  We  make  the  assumption 
that  the  start  and  stop  of  the  earthquake  are  all  included  within  the  initial  P-wave  and  the  source  s 
radiation  is  described  by  Savage's10  model.  If  this  is  so.  and  since  the  event  discriminates 
(M  =5  5  M  for  Longshot  =  4.0),  we  cannot  invoke  either  source  dimensions  or  an  elongated 
source-time  function  as  leading  to  discrimination.  We  could,  of  course,  have  a  radically  differ¬ 
ent  source-time  function  between  the  two  events,  such  as  a  step  for  an  earthquake  and  an  impulse 
for  an  explosion,  but  nowhere  in  the  literature  is  there  any  good  evidence  that  explos.ons  are  dom¬ 
inantly  impulsive  in  displacement.  Indeed,  it  appears  that  when  many  explosions  have  been  studied 
as  a  suite  11  it  has  been  impossible  to  make  a  clear-cut  statement  about  source-time  functions,  as 
some  appear  to  have  small  impulsive  terms  while  others  appear  to  be  simple  steps.  (F.lson  and 
Lande  reached  a  similar  conclusion  in  the  SATS  dated  31  December  1970,  DDC  AD-718971.)  All 
discriminate  easily,  however.  1  conclude,  therefore,  that  it  would  be  difficult  to  sustain  an  argu- 
rne.  1  that  source-time  functions  or  source  dimensions  control  discrimination. 

D.  Davies 


E.  THE  R  PHASE  FROM  EXPLOSIONS  AT  NEAR  DISTANCES 

An  observational  study  of  the  R  phase  has  been  made  using  two  stations  within  300  km  of 
explosive  sources  at  the  Nevada  Test  Site.  R  is  taken  to  represent  a  high-frequency  Rayleigh- 
type  motion  which  propagates  in  crustal  wave  guides.  The  chief  purpose  of  the  experiment  was 
to  determine  if  R  could  be  used  to  discriminate  between  contained  underground  explosions  and 
those  used  for  excavation  purposes.  The  maximum  amplitude  of  this  phase  was  measured  on 
short-period  vertical  instruments  at  Mina.  Nevada  (MN-NV,  A  =  240km)  and  Kanab,  Utah 
(KN-UT.  A  =  290km)  from  explosions  detonated  under  various  conditions.  This  maximum  am¬ 
plitude  is  usually  recorded  35  to  50  sec  following  the  initial  P-wave  arrival  at  these  distances, 
and  is  not  a  well-defined  phase  but  consists  of  a  series  of  high-frequency  excursions  of  larger 
amplitude  than  the  preceding  motion.  The  amplitudes  (peak  to  trough)  of  the  initial  (first  full 
cycle)  P  arrival  (P  )  and  the  long-period  fundamental  mode  Rayleigh-wave  motion  <LR)  were 
also  measured  for  comparison  with  the  R  measurements.  The  amplitude  measurements  were 
corrected  for  instrument  response  and  the  average  using  the  two  stations  computed  directly. 

Only  explosions  with  source  conditions  listed  by  Springer  and  Kinnaman  were  used  in  this  study. 
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tn  Fig.  1-12,  the  average  surfa  e  wave  amplitudes  1,^  are  potted  against  I’n  for  some 
24  explosions;  various  symbols  are  used  to  denote  the  source  medium  and  location.  The  ex¬ 
plosions  termed  "Plowshare  excavation"  by  Springer  and  Kinnanian12  are  circled;  the  explosion 
described  as  a  "Plowshare  experiment"  is  enclosed  in  a  square.  Given  the  scatter  of  the  data 
in  Fig.  1-12,  there  is  nothing  remarkable  to  define  the  excavations  except  that  their  surface  wave 
amplitudes  may  be  larger,  relative  to  P  ,  than  some  of  the  contained  explosions.  In  Fig.  I- 1  5, 
we  ha,  ted  II  r  vs  Mere,  11  for  at  least  one  of  the  excavations  (named  Sedan,  detonated 

in  alluvium  at  Yucca  Flat)  is  significantly  lower  relative  to  I,R;  however,  the  same  is  true  for 
the  explosions  named  Palanquin  and  Scroll,  not  specifically  described  as  excavations.  H  for 
the  other  excavation  event  Schooner  is  low  hut  not  greatly  so.  The  most  interesting  aspect  of 
Fig.  1-13  is  the  apparent  difference  in  dependence  of  upon  I.jj  when  compared  with  Fig.  1-12 
where  the  data  may  be  roughly  approximated  by  a  line  with  a  slope  of  0.7  while  a  less  steep  line 
w(ju.  .  ,  re  sc  »  the  data  of  Fig.  1-1  3.  If  we  take  to  be  a  measurement  of  relative  yield,  we 

must  conclude  mat  explosions  of  larger  yield  do  not  generate  as  much  H  relative  to  P  as  do 
smaller  ones,  clue  either  to  their  size  or  their  greater  depth. 

In  Fig.  1-14,  we  test  the  latter  hypothesis  hy  plotting  the  ratio  y  (where  y  -  vs 

source  depth.  It  is  clear  from  this  figure  that  y  increases  with  depth,  the  deeper  events  being 
relatively  more  efficient  in  coupling  energy  into  1.^  than  into  11^.  Additionally,  the  excavation 
explosions  are  now  clearly  separated  fi  mi  the  contained  explosions  at  the  same  depth;  however, 
Palanquin  and  Scroll  have  y  values  nearer  the  excavation  explosions  than  to  the  contained  events. 
We  cannot  assume  Palanquin  was  an  excavation  device,  being  described  only  as  a  "Plowshare 
experiment"  in  the  source  list.  The  source  medium  of  Scroll  is  described  as  "tuff  (volcanic  ash)," 
which  may  account  for  its  high  value  of  T. 

Of  course,  in  a  relevant  application  of  those  empirical  results,  the  depth  of  the  explosion 
will  be  unknown.  In  Fig.  1-15  we  have  plotted  7  vs  P  ,  both  quantities  which  are  measureable 
from  a  seismogram.  Here  Schooner,  the  Pahute  Mesa  excavation  in  tuff,  plots  only  slightly 
higher  than  explosions  contained  in  tuff  at  Yucca  Flat;  however,  the  value  of  y  for  Sedan  is  about 
an  order-of-magnitude  greater  than  a  contained  alluvium  explosion  with  the  same  P  value.  Again, 
y  for  Scroll  is  high  and  in  this  case  that  of  Palanquin  is  not  anomalous. 

To  summarize  and  conclude  briefly,  we  have  measured  the  amplitudes  of  Pn<  H^,  and  L1I 
at  two  stations  within  300  km  of  24  explosions.  The  ratio  i,fj/l?g  appears  to  increase  with  depth, 
and  excavation  explosions  show  a  value  of  this  ratio  greater  than  those  of  contained  explosions 
at  a  similar  depth.  This  ratio,  relative  to  Pn  amplitude,  is  markedly  higher  for  an  excavation 
explosion  in  alluvium  compared  with  contained  explosions  at  the  same  site.  The  excavation  ex¬ 
plosion  in  tuff  does  not  show  this  feature  to  such  an  extent;  however,  no  data  for  a  contained  tuff 

event  of  the  same  P  amplitude  have  hern  considered  thus  far. 

n  J.  Filson 

L.  Lande 
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CENTERED  *T  44  O*  N  77.»*E  (KA) 


Fig.  1-1.  Mop  of  Eurasia  showing  locations  of  12  WWSSN  stations  (squares)  employed 
to  anclyze  7  events  (crosses)  studied  in  this  SATS.  Refer  fa  Table  1-1  for  event  data. 
Great  circle  paths  between  events  and  stations  are  approximately  straight  lines. 
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TB  SHL  St  5*  4  9.33 


TB  CMO  3k  \f  49:3  3 


Fig.  !-"*  Vertical  components  of  Rayleigh  waves  from  events  KA,  TB,  KB,  and  KD  at  comparable 
distances  and  gains.  KD's  seismogram  has  amplitudes  that  are  slightly  smaller  than  average  con¬ 
tinental  Eurasian  earthquake  with  same  m(j.  High-frequency  Rayleigh  wave  from  event  KA  rides 

atop  interfering  low-frequency  Rayleigh  wave  from  event  in  Tonga  (m,  =6.0). 
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EVENT  KA  5*1-69  AT  KBL  6h 


Fig.  1-4.  Three  components  of  surface  waves  generated  by  event  KA  os  recorded  by  Kabul. 
Energy  arriving  with  group  velocity  of  3.34  km/sec  may  be  either  Lr-ve  or  nigher-mode 
Rayleigh  waves. 
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Fig.  1-5.  Short-period  signals  from  event  KA  as  recorded  by  arrays  at  LASA,  Yellowknife, 
and  Warramunga,  showing  appaient  depth  phase  and  waveform  simplicity. 


Fig.  1-6,  (a)  l(t),  the  impulse  response  of  long-period  seismometer,  (b)  I  (t),  the  inverse 
impulse  response  computed  using  Eq.(l-2).  (c)  Result  of  M(t)  *  l(t).  If  computation 
of  l”l(t)  were  exact,  this  should  be  a  delta  function. 
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tin:  i  hi  i 


Fig.  1-7.  (o)  A  long-period  P-wove  os  recorded,  S(t).  (b)  Ground  motion  P(t)  computed 
using  Eq.(l-3),  i.e.,  I  (t)  *  5(t).  (c)  Ground  motion  P(t)  computed  by  direct  inversion 
of  Eq.  (1-1). 
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Fig.  1-9.  Azimuthal  variation  of  deconvolved,  vertical  long-period  P-woves  with  respect  to  fault  pli 
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BKS 

Fig.  1-10.  Deconvolved  P-v*aves  at  ATU 
and  BKS  showina  e'ongation  at  BKS  taken 
to  be  result  of  fault  propagation. 
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Fig.  1-11.  Short-period  LASA  P-wave  signals  from  seismic  events 
in  Aleutians;  (a)  and  (b)  traces  are  P  and  PcPfrom  Komandorsky 
event  of  20  January  1969,  and  (c)  trace  is  PcP  from  Longshot. 
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Fig.  1-12.  Average  of  F  amplitude  at  MN-NV  and  KN-UT  plotted  vs  average  L  amplitude. 
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Fig.  1-13.  Averoge  maximum  R  amplitude  plotted  vs  overage  amplitude. 
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K  ' 1  ■  "fit 

Fig.  1-14.  Average  value  of  J,  where  7  =  L^/R  ,  plotted  vs  depth  of  explosion. 
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Fig.  1-1.*;.  y  plotted  vs  P  . 
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II.  PROPAGATION 


A.  THF  EFFECT  OF  LATERAL  HETEROGENEITY 
UPON  SURFACE  WAVE  PROPAGATION 


In  several  recent  SATS,  the  problem  of  body  wave  propagation  in  laterally  heterogeneous 
earth  models  has  been  considered,  with  a  view  to  predicting  how  variations  in  the  eartii  may 
distort  seismic  observations  and  complicate  the  di  icriminat ion  problem.  We  found  that  pro¬ 
nounced  focusing  and  defocusing  (or  shadowing)  of  energy  can  occur  and  will  introduce  systematic 
errors  into  measurements  of  body  wave  magnitude.  Our  knowledge  of  the  earth's  internal  struc¬ 
ture  is,  however,  not  yet  good  enough  to  enable  us  to  make  quantitative  predictions  of  these  ef¬ 
fects  in  general,  although  we  can  probably  predict  shadow  locations  for  sources  in  island  arcs 
with  reasonable  accuracy 

For  surface  waves,  similar  propagation  effects  occur,  and  they  are  more  amenable  to  nu¬ 
merical  calculation  and  prediction,  since  surface  wave  velocities  are  functions  of  two  spatial 
coordinates,  rather  than  tiirce,  and  are  correlated  with  known  geologic  and  tectonic  features. 
Furthermore,  our  seismic  discrimination  capability  is  currently  limited  by  the  surface  wave 
detection  threshold,  so  reliable  path  corrections  would  be  of  great  value.  Tire  ability  to  reliably 
calculate  probation  effects  might,  in  fact,  lower  the  surface  wave  detection  threshold  by 
facilitating  the  design  of  matched  filters. 

To  analyze  surface  wave  propagation,  we  have  used  geometrical  ray  theory  -  much  as  is 
commonly  done  for  body  waves.  For  surface  waves,  howiver,  some  additional  complications 
arise  because  of  dispersion.  In  particular,  in  a  iieterogeu  us  dispersive  medium  the  group 
and  phase  velocities  do  not  generally  have  the  same  directh  n.  lly  considering  the  Interference 
of  two  plane  waves  with  slightly  different  frequencies  w  ami  wavenumber  vectors  k,  it  can  be 
seen  that  the  group  slowness  is  given  by 


Since  k  is  a  function  of  the  initial  parameters  of  a  ray  and  w,  ai  well  as  the  observation  point, 
the  calculation  of  eitiier  the  magnitude  or  the  direction  of  the  group  slowness  is  not  easy.  How¬ 
ever,  since  for  any  given  source 

k(r,  u)  Mr,  u)  ii(F,  u) 


where  ii  is  a  unit  vector  normal  to  the  (dta.se  wavefronts,  and  ?  is  the  position  vector,  we  have 


d£ 

Du 


Dk  -  .an 

t;  -  n  i  k  - — 
<Ju  o  tv 


Since  n  is  a  unit  vector  (Dn/Ou)  •  n  0  and 

Dk  -  Dk 
Du  n  Du 


This  result  says  that  although  tile  actual  group  dowries*  vector  is  difficult  to  calc  date,  its  com¬ 
ponent  in  the  ray  direction  is  the  conventional  (and  easily  calculable)  group  slov  ness  for  a  homo¬ 
geneous  medium.  The  group  delay  could  be  obtained  by  integrating  the  group  si-  wiles*  along 
any  pith,  but  for  die  phase  ray  path,  this  task  is  greatly  simplified. 


Preceding  page  blank 
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One  of  the  most  striking  pieces  of  evidence  that  lateral  heterogeneity  tsignifican  ly  affects 
surface  wave  propagation  is  the  discovery  by  Capon1  that  Rayleigh  waves  observed  at  I  ASA  often 
avel  in  directions  greatly  different  from  the  great  circle  azimuth  to  the  event  that  generated 
t  torn.  I'igu  e  11- t,  for  example,  shows  the  paths  inferred  by  Capon  to  have  been  followed  by 
the  Itayloigh  waves  from  an  earthquake  in  the  Kuril  Islands.  Capon  felt  that  most  of  the  observa¬ 
tions  he  made  could  be  explained  by  refraction  and/or  reflection  at  continental  margins.  In  order 
to  check  this  assertion,  a  regionalized  earth  model  consisting  only  of  continents  and  oceans  was 
constructed  and  used  for  ray-tracing  calculations.  The  phase  velocities  for  oceans  and  continents 
were  assumed  to  be  1.50  and  3. OK  km/ sec  at  20-sec  periods,  and  the  model  was  specified  .n  terms 
ot  numerical  values  on  a  S'  X  0  latitude-longitude  grid.  The  resulting  ray  paths,  shown  in 
fig.  11-2,  look  remarkably  similar  to  those  inferred  by  Capon,  and  strikingly  confirm  the  imprr- 
tanco  of  continental  margins  which,  in  diis  case,  have  deflected  the  rays  by  total  internal  reflec¬ 
tion.  (  nrrently,  work  is  being  done  to  increase  the  amount  of  detail  in  the  model,  to  see  how 
well  propagation  effects  can  be  predicted  quantitatively. 

H.  R.  .lulian 

B.  TRAVEL-TIME  STUDY  USING  DEEP  EARTHQUAKES 

AND  LATERAL  HETEROGENEITIES  IN  THE  DEEP  MANTLE 

An  attempt  has  been  made  to  refine  the  travel-time  curve  for  telcscismic  l’-waves  us  ng 
data  from  deep-focus  earthquakes  This  approach  was  used  in  order  lo  minimize  the  bias  intro¬ 
duced  into  travel  limes  by  lateral  heterogeneities  in  the  upper  mantle,  which  are  |uriicularlv 
severe  neai  island  arcs,  where  seismic  activity  is  concentrated.  The  raw  data  used  (nearly 
•1  '(10  arrival  times  from  earthquakes)  were  taken  from  bulletins  of  the  International  Sels- 
iiiohigtt.il  C  intro  and  the  National  Ocean  Snrvev  (NOS),  and  were  winnowed  to  eliminate  gross 
errors  by  comparing  data  for  several  earthquakes  in  each  sotirre  region.  An  iterative  procedure 
of  relocating  the  events,  modifying  the  travel -time  curve  and  station  correct  ini  s,  relocating 
the  events,  etc.  was  then  used  to  determine  finally  the  average  travel-time  curve  and  a  set  of 
station  corrections.  The  method  will  be  described  In  detail  in  a  paper  no\  t.eing  prepared. 

El  cure  11- 1  shows  the  travel -time  curve  derived  for  a  focal  depth  of  5*>  *m.  An  ambiguity 
still  remains  with  regard  to  tie  absolute  values  of  the  times,  since  the  origin  times  of  the  <artli- 
quakes  cannot  be  determined  inde|smdemlv.  We  lave  used  data  from  Nevada  Test  Site  explosions, 
where  the  stria  tnre  of  the  upper  mantle  is  fairly  well  known,  to  estimate  this  *!).('.  coni|>onciii,” 
and  found  that  If*  see  should  be  subtracted  from  the  values  In  I'lg.  11-3.  This  Is  the  value  which 
would  applv  to  a  deep-focus  earthquake  ill  Nevada. 

Hie  general  shape  of  our  curve  Is  similar  to  those  found  by  other  workers  (e.g.,  Herrin. 

Si!  jJLTt  although  our  absolute  times  are  nl-out  I  tn  I  5  sei  later  than  llerrln  s.  Ilevond  *5  ,  how¬ 
ever.  tit*’  trend  of  our  <  ttrve  is  decide  11  v  different  from  mwt  others  ami,  moreover,  is  subje<  t 
o  lane  regional  varl.itp.ns.  the  tin  reused  scatter  of  the  data  hevond  H<*  is  evident  In  I  Ig  11-3. 
We  suggest  that  this  effect  Is  i  a  used  In  regional  variations  hi  tin’  strut  litre  of  the  deep  mantli 
•Statistical  tests  establish  clejrlv  tb.it  the  variations  are  .systematic,  not  ran  bun,  and  tin*  only 
systematic  error  which  might  be  stronglv  distance  dependent  t.s  one  caused  bv  variations  near 
the  liottomtiig  point*  of  the  rav  1  tgurn  11- »  shows  region*  lo-Ueved  to  have  |,|C|,  and  low  ve|..«  - 
Hie*  111  die  deep  man'.lu.  Miitil  ii  analysis  was  performed  for  other  depth*  and  vet  v  few  .,j  tnii- 
cant  variations  could  be  f  mu  !  ttoik  1*  curr**iitl,  being  >uri  (o  J  out  t  >  refine  tlc-4«  rftMtam  , 
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to  obtain  coverage  of  more  of  the  earth,  and  to  include  other  properties  such  as  shear  veloeity 


and  attenuation. 


B.  It.  Julian 
M.  K,  Sen  Gupta 


C.  THE  ARRAY  LOCATION  PROBLEM 

The  standard  technique  for  locating  with  an  array  such  as  LASA  is  to  fit  a  least-squares 
plane  wave  to  the  arrival  times  at  many  sensors  and  hence  to  determine  dT/dA  and  azimuth  for 
the  incoming  wave.  We  soon  realized  that  this  approach  was  unsatisfactory  in  that  the  signal 
arrival  times  often  did  not  lie  close  to  the  plane  predicted  for  earthquakes  from  known  locations. 
Two  separate  components  may  be  identified  for  this  misfit.  An  incoming  plane  wave  is  unlikely 
to  traverse  the  crust  beneath  LASA  without  some  degree  of  corrugation  being  acquired  from 
small-scale  geological  heterogeneities.  On  the  other  hand,  the  best-fitting  least-squares  plane 
wave,  in  whieh  we  attempt  to  average  out  these  corrugations,  may  also  have  a  dip  (or  dT/dA) 
and  azimuth  different  from  that  predieted  by  standard  travel-time  tables,  It  is  this  effeet  which 
we  consider  here.  Of  course,  routine  operation  of  the  array  for  bulletin  production  purposes 
circumvents  the  scientific  problem,  sinee  empirieal  corrections  may  be  applied.  However,  when 
an  event  in  a  new  area  has  to  be  examined,  it  may  be  impossible  to  interpolate  corrections  ade¬ 
quately.  Thus,  we  need  to  have  a  good  understanding  of  the  cause  of  array  mislocation. 

Three  eauses  may  lie  considered: 

(1)  If  there  is  broad-scale  structure  beneath  I  .ASA,  such  as  a  dip  of  the 
Molio,  there  will  be  a  smooth  bias  in  the  mislocation.  This  has  been 
extensively  studied,  notably  by  Greenfield  and  Sheppard 3  and  by 
Fngdaltl  and  Felix.'* 

(2)  If  the  reference  dT/dA  tables  (Jeffreys-Bullen)  were  incorrect  at 
some  distances,  the  clip  (but  not  the  azimuth)  of  the  planes  would  lie 
affected  for  all  events  at  these  distances. 

(B  If  there  is  lateral  heterogeneity  in  the  earth,  rays  will  not  necessarily 
trace  great  circle  paths,  and  dT/dA  and  azimuth  anomalies  will  be 
equally  likely.  The  scale  of  these  anomalies  may  be  inferred  from 
LASA  data.  Near-sotiree  structure  falls  within  this  category. 

In  order  to  attempt  to  isolate  these  three  effects,  we  show  in  Fig,  11-5  an  "array  diagram." 
Tile  radial  coordinate  is  dT/dA,  (lie  azimuthal  eoordinate  is  azimuth.  Bach  event  studied  is 
represent  'd  by  an  arrow,  the  head  of  which  indicates  the  expected  position  on  the  diagram  if  the 
NOS  location  were  correct,  if  the  .1-1$  tables  for  dT/dA  were  accurate,  if  there  were  no  lateral 
heterogeneity,  and  if  there  were  no  varying  structure  under  the  array.  The  tail  of  each  arrow 
represents  the  measured  dT/dA  and  azimuth  (hoping  that  corrugations  on  the  wave  front  will 
introduce  no  unusual  Idas).  Thus,  the  arrow  length  and  direction  indicate  the  extent  to  which 
the  conditional  clauses  above  arc  satisfied. 

The  arrav  diagram  (thinned  somewhat)  is  for  a  large  population  of  events  recorded  at  LASA 
from  I  If, 7  through  IV/*.  \\  e  emphasize  that  no  prior  corrections  have  been  made  in  thedT/dA 
an  azimuth  measurements  these  arc  raw  data.  The  problem  is  to  isolate  the  different  causes 
for  arrows  not  to  be  of  zero  length.  We  shall  assume  at  the  outset  that  the  NOS  location  is 
correct  winch  is  probably  an  acceptable  assumption  at  the  precision  to  which  we  are  working. 
The  three  remaining  sources  of  bias  then  affect  the  array  diagram  in  somewhat  different  ways. 
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A  smooth  distribution  of  arrows,  all  approximately  of  equal  length  and  all  pointing  in  the 
same  direction  would  characterize  bias  introduced  beneath  the  array.  A  distribution  of  arrows 
i  which  radial  terms  dominated  and  were  of  equal  magnitude  at  all  azimuths  would  characterize 
ias  arising  from  use  of  an  incorrect  (globally)  dT/dd  curve.  A  distribution  of  greater  complexity 
would  point  to  near-source  or  propagation  path  bias,  including  equally  azimuthal  and  dT^A 
anomalies. 

It  is  clear  from  Fig.  11-5  that  the  situation  is  immensely  complex.  Can  we  extract  even  a 
simple  sti  ucture  beneath  the  array  from  this?  We  show  in  Fig.  11-6  the  distribution  of  arrow 
lengths  and  directions,  and  to  the  extent  that  this  is  a  random  sampling  we  may  conclude  that 
there  is  a  mean  bias  to  the  north  of  0.25  sec/deg.  We  tentatively  associate  this  bias  with  the 
smooth  effects  that  would  be  introduced  by  structure  beneath  the  array  and  which  would  move 
the  tail  of  each  arrow  0.25  unit  south. 

T'ic  arguments  for  the  relative  merits  of  heterogeneity  near  source  and  along  propagation 
path  are  complicated  and  will  be  listed  at  length  in  a  forthcoming  publication.  The  array  struc¬ 
ture  bias  is  smooth  and  we  are  discussing  rapid  fluctuations  in  the  array  diagram.  Thus,  we 
are  confident  that  our  assertions  are  truly  about  the  source  and  propagation  path.  We  see  no 
compelling  evidence  that  there  is  a  global  dT/dA  anomaly  shown  by  the  array.  As  a  result  of 
detailed  analysis  of  Fig.  11-5,  we  conclude: 

(1)  There  is  lateral  heterogeneity  in  the  earth's  mantle  with  good 
indications  that  near  the  core  mantle  boundary  there  may  be 
horizontal  structures  of  100-km  extent  or  less. 

(2)  The  amplitude  of  these  anomalies  suggests  velocity  fluctuations 
of  perhaps  5  percent,  but  a  full  interpretation  will  be  very 
difficult  with  a  data  set  restricted  by  seismicity. 

(“>)  Regions  of  marked  lateral  variation  near  the  core  boundary 

include  areas  beneath  Hawaii,  the  Galapagos  and  the  mid-Atlantic 
Ridge  north  of  Iceland.  This  result  may  be  of  use  in  examining 
Morgan's^  hot-spot  hypothesis. 

(4)  Lateral  heterogeneity  is  by  no  means  confined  to  the  deepest 
mantle.  There  is  some  evidence  that  sub-continent  and  sub-ocean 
structural  differences  may  be  detectable, 

(5)  Many  anomalies  remain  to  be  understood,  particularly  in  the 
north-west  quadrant  as  Chinnery^  has  described. 

(6)  In  any  study  using  arrays  to  examine  the  velocity  structure  of 
the  earth,  it  is  imperative  to  see  both  dT/dA  and  azimuth  of 
incoming  signals,  as  azimuth  deviations  can  be  as  great  as  10° 
and  bear  as  much  useful  information  as  dT/dA  anomalies  —  indeed, 
the  two  cannot  be  separated  in  a  heterogeneous  earth. 

From  an  operational  point  of  view,  we  remark  that  the  location  of  an  event  in  a  new  and 
aseismic  area  will  present  problems  not  necessarily  alleviated  by  extrapolation  of  existing 
results.  We  are  n>-.v  starting  to  examine  NORSAR  in  a  similar  way. 


I).  Davies 
R.  M.  Sheppard 
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Fig.  11-1.  Rayleigh-wove  propagation  paths  for  earthquakes  in  Kuril  Islands 
inferred  by  Capon'  from  observations  made  using  LASA. 
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Fig.  11-3  P-wave  travel-time  curve  derived  in  this 
study,  expressed  in  terms  of  syslemafic  deviation 
from  values  in  Jeffreys-Bullen  tab  es  (focal  depth 
550  km).  Mean  values  and  standard  errais  af  mean 
(vertical  bars)  in  2°  intervals  are  shown,  as  are 
smoothed  values  at  5°  intervals  (squares). 


DISTANCE  (d*5 


Fig .  II -4,  Su  ested  regions  of  low  (L)  and  high  (H)  campressional  velocities  in  deep  mantle  (2600  km  depth 
<2900  km).  ry  paths  used  ,hinned)  are  shown. 
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Fig.  11-5.  Array  diagram  far  LASA.  Outer  ring  has  radius  af  10  sec/deg.  Core  data  are  not  considered 
in  this  study. 
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We  now  wish  to  define  the  basic  quantities,  derived  from  the  Hit  analysis,  which  are  impor¬ 
tant  in  the  interpretation  of  the  P-wave  anomaly  data.  These  quantities  are  the  azimuth  devia¬ 
tion  (AZl)E'  ),  the  slowness  correction  (SLOCOIt),  and  the  relative  power  expressed  in  decibels 
(PDBR).  The  definition  for  AZDEV  is 


AZ DEV  =  MEASURED  AZIMUTH  -  TRUE  AZIMUTH 

Tiie  ti  ue  azimuth  is  assumed  to  be  that  obtained  from  the  CCS  location  for  the  epicenter  of  the 
event.  In  -  uefinition  for  SLOCOR  is  taken  to  be  consistent  with  that  given  by  Niazif*  The  def¬ 
inition  for  SLOCOR  is 


SLOCOR 


SLOWNESS  IN  MANTLE 
SLOWNESS  AT  SURFACE 


where  the  slowness  at  the  surtaee  is  the  reciprocal  of  the  horizontal  phase  velocity  obtained  as 
discussed  previously  from  the  HR  analysis,  and 

SLOWNESS  IN  MANTLE  =  sin  (L.)/8.2 

IV 1 


where  8.2km/see  is  the  assumed  coi  .pressional  velocity  in  the  mantle,  and  L,  is  the  angle  of 
incidence  of  the  P-wave  in  the  mantle*  The  angle  is  obtained  from  Appendix  V  of  lliehter, 
using  the  appropriate  epicenter -to- 1  .ASA  distance  for  the  event  and  noting  that  a  crust-to-mantle 
conversion  must  be  made  of  the  angles  tabulated  by  Richter. 

The  definition  for  PDBR  is 


PDDR  =  10  log1Q  (PMAX/POWER) 

where  POWER  is  the  power  (in  m g2/H^t  measured  using  the  HR  method,  and  PMAX  is  the  max- 

2 

imum  power  measured  (in  mg  /Hz)  among  the  9  subarrays,  for  each  event.  It  should  lie  stressed 
that  there  will  be  279  independent  measurements  for  AZDEV,  SLOCOR,  and  PDBR.  That  is, 
there  are  ,1  events  for  which  these  quantities  have  been  measured  within  the  9  subarrays  com¬ 
prising  the  A-,  B-,  C-rings  at  1  .ASA, 

A  graphic  illustration  of  the  data  for  AZDEV  and  SLOCOR  is  given  in  Pig.  111-2  for  the 
26  February  1967  Eastern  Kazakh  event.  In  this  figure,  the  diameters  and  locations  of  the  circles 
represent  the  apertures  and  locations  of  the  various  subarrays,  and  the  diameters  of  the  circles 
are  drawn  along  the  true  azimuth  of  the  event.  The  arrows  point  in  the  direction  of  propag.  lion 
of  the  wave,  and  the  length  of  the  arrow  is  proportional  to  SLOCOR,  at  the  various  subarrays. 

The  strikingly  large  changes  in  AZDEV  and  SLOCOR  between  subarrays  which  are  located  within 
about  one  wavelength  of  the  P-wave  from  each  other  should  be  noted.  This  fact  is  extremely 
important  in  the  interpretation  of  the  data  in  terms  of  a  possible  structure  within  the  earth  which 
could  produce  such  anomalies. 

It  is  instructive  to  plot  the  data  for  AZDEV,  SLOCOR,  and  PDBR  as  a  function  of  the  true 
azimuth  ot  the  epicenter  of  the  event,  as  depicted  in  Kigs.  HI  —  3(a)  through  (c),  respectively.  In 
these  fij  ures,  an  average  curve  is  drawn  which  represents  the  average  of  the  data,  for  the  9  sub¬ 
arrays,  at  each  azimuth  or,  equivalently,  for  each  event.  It  should  be  noted  that  SLOCOR  and 
PDBR  have  been  limited  to  be  less  than  3.0  and  15.0dB,  respectively.  The  relatively  large 
values  for  AZDEV  and  SLOCOR  that  occur  in  Pigs.  II I- 3(a)  and  (b),  respectively,  are  for  the 
core  phases  which  were  used  in  the  analysis.  This  is  to  be  expected,  since  the  angle  of  incidence 
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for  those  phases  is  small,  so  that  the  AZDEV  and  Sl.OCOR  for  them  should  be  large,  cf,  Niazi 
An  interesting  observation  to  make  at  this  point  is  that  the  average  curves  for  AZDEV  and  SLOCOR 
m  Figs.  Ill- 3(a)  and  (b),  respectively,  appear  to  vary  about  zero  and  unity,  respectively,  with  a 
certain  amount  of  regularity.  For  example,  in  Fig.  !ll-3(a)  the  average  curve  for  AZDEV  appears 
to  be  negative  for  azimuths  between  150°  and  330°  and  positive  otherwise;  that  is,  the  average 
curve  in  the  figure  makes  a  transition  between  negative  and  positive  values  at  an  azimuth  of  330°. 

In  addition,  the  average  curve  for  Sl.OCOR  in  Fig.  Ill- 3(b)  appears  to  be  less  than  unity  for  az¬ 
imuths  between  60°  and  240°  and  is  greater  than  unity  otherwise.  According  to  the  results  of 
Niazi,"  this  type  of  behavior  is  highly  suggestive  of  the  effect  introduced  by  a  plane-dipping  inter¬ 
face  across  which  there  is  a  change  in  eompressional  velocity. 

We  attempted  to  find  the  dip  and  direction  of  strike  of  a  plane-dipping  interface  that  would 
provide  a  best  least-squares  fit  to  the  AZDEV  and  Sl.OCOR  data.  A  velocity  ratio  of  0.75  was  as¬ 
sumed  and  the  procedure  given  by  Niazi^  was  employed  to  compute  the  AZDEV  and  Sl.OCOR,  for 
an  assumed  plane  interface,  for  each  event.  We  found  that  the  dip  angle  of  such  an  interface  was 
11°  and  that  the  azimuth  of  the  strike  direction  was  60°,  This  agrees  quite  well  with  the  average 
dip  angle  and  direction  of  strike  of  the  corrugated  surface  proposed  by  Greenfield  and  Sheppard1 
to  explain  the  P-wave  time  anomaly  data  at  I  ASA.  Additional  work  has  been  done  to  determine 
a  structure  within  the  earth  to  explain  the  P-wave  anomaly  data,  and  the  results  will  be  presented 

U 

in  a  forthcoming  publication,  A  brief  discussion  of  some  of  these  results  will  now  be  presented. 

A  model  has  been  proposed  to  explain  the  cause  of  the  P-wave  anomaly  data,  in  which  the 
crust  and  upper  mantle  structure  under  I  ASA  can  be  characterized  as  a  random,  or  Chernov- 

9 

type,  medium.  A  statistical  analysis  of  the  data  for  the  phase  and  amplitude  fluctuations  of  the 

n 

P-wave  has  been  made  which,  along  with  some  theoretical  results  due  to  Chernov,  was  used  to 
determine  the  properties  of  the  random  medium.  We  found  that  the  correlation  distance  of  the 
random  medium  is  12  km  and  that  the  standard  deviation  of  the  index  of  refraction  is  1 . 9  percent. 
These  random  inhomogeneities  exist  in  the  crust  and  upper  mantle  of  the  earth  down  to  a  depth 
of  about  1  36  km. 

It  is  important  to  point  out  that  a  single,  or  multiple,  corrugated  interface  cannot  be  used 

to  explain  the  source  of  the  P-wave  anomalies  because  such  a  structure  would  lead,  according 

to  Chernov's  theory,  to  a  high  correlation  between  the  amplitude  and  phase  fluctuations,  and  to 

a  value  of  amplitude-relative-to-phase  fluctuation  power  ratio  close  to  zero.  These  conditions 

are  not  found  to  be  true  in  the  statistical  analysis  of  the  data.  Hence,  a  model  involving  a  single, 

or  multiple,  corrugated  interface  is  inadequate  to  explain  the  observations.  It  should  be  men- 

10 

tioned  that  a  similar  approach  to  this  problem  has  been  considered  recently  by  Aki. 

.1,  Capon 


B.  CRUSTAL  SCATTERED  S-WAVES  UNDER  LASA 

The  nature  of  the  laterally  heterogeneous  structure  of  the  crust  and  upper  mantle  beneath 
[ASA  that  causes  the  first-order  time,  amplitude,  and  waveform  anomalies  associated  with 
short -period  tob. seismic  P -waves  is  being  investigated  by  observing  the  scattered  S-waves  from 
emerging  l’-waves.  Such  a  study  is  presently  confined  to  data  obtained  from  the  3-componeut, 
short-period,  digitally  recorded  seismograph  at  subarray  1)2  (see  Ref.  11).  |The  data  are  re¬ 
corded  on  all  standard  high-rate  format  1A.SA  tapes  In  record  words  19  (Vertical),  40  (North), 
and  61  (East)  after  26  dtine  1968.] 
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12-1-1 

Previous  investigators  have  used  only  the  short-period  P-wave  motion  recorded  on  the 

vertical  instruments,  lu  obtain  a  model  based  on  such  data,  one  must  know  the  first-order 
P-wave  velocity  structure  beneath  the  array.  Observing  the  total  motion  has  several  distinct 
advantages.  I'ht  S-veloeity  as  a  function  of  depth  is  not  as  well  known  as  the  P-velocity  function; 
however,  the  location  of  a  source  of  S  will  lie  more  accurate  when  determined  by  the  first-order 
^-P  tithe  measurement  than  when  derived  by  travel-time  residuals  alone.  Further,  since  S  and 
P  '  -Jentiallv  normally  polarized,  a  simple  coordinate  rotation  separates  them.  The  rela¬ 
tive  amplitude  am!  frequency  content  of  the  two  waves  provide  information  on  the  impedance 
contrast  and  sharpness  of  any  boundary  acting  as  an  apparent  source.  The  same  functions  de¬ 
termined  over  various  distances  ami  azimuths  determine  the  strike  and  dip  of  the  converting 
boundary.  Moreover,  the  shape  of  a  converted  S-wave  will  essential!-  he  the  same  as  that  of 
the  exciting  P-wave,  thus  the  observed  differences  between  tin?  S-  and  P-wavcl'orms  are  cun- 
si.^  i-d  in  be  the  result  of  a  process  occurring  between  the  apparent  source  and  the  receiver, 
Cuiisequen*  I v,  dull  part  ol  the  seismogram  due  to  crustal  reverberation  should  be  separable  from 
the  source  component, 

Data  from  a  typical  event  are  shown  in  Fig.  Ill— I.  The  raw  data  in  the  first  three  traces 
show  P  and  pP  from  an  event  located  in  the  North  Atlantic  at  a  distance  of  f>8"  and  an  azimuth 
ol  b2  ",  The  bottom  two  traces  show  the  rotated  components  containing  SV  and  SI!  motion.  The 
P-oomponent  is  approximately  the  same  as  the  V -component  At  a  lag  of  1  to  2  see  behind  P 
and  pP  is  an  S-wave  arrival.  Such  an  oflset  constrains  the  apparent  source  of  S  to  lie  at  a  depth 
of  about  3km,  the  boundary  between  the  sediments  and  basement  at  i ASA. 15  Tin-  fact  that  the 
S-wave  from  pP  has  not  been  completely  polarized  at  the  same  angle  as  the  S-wave  from  P  in¬ 
dicates  that  P  and  pP  are  emerging  from  different  azimuths.  Relative  particle  motion  on  the  SV 
and  SI  1  components  of  the  S-wave  due  to  pP  determine  that  the  apparent  azimuth  of  pP  is  about 
20° N  of  the  apparent  azimuth  of  P.  The  frequency  content  and  amplitude  of  the  S-waves  (0.35  of 
the  P-wave  amplitude)  are  consistent  with  a  tectonic-sedimentary  boundary  that  is  planar,  sharp, 
and  has  a  high  impedance  contrast  (refraction  studies*  show  a  velocity  contrast  of  about  2  for  this 
boundary;  however,  no  conclusion  on  the  sharpness  or  uniformity  of  the  boundary  can  lie  obtained 
from  refraction  data).  The  strike  and  dip  of  this  boundary  can  be  deiermined  by  observing  the 
S-waves  it  scatters  for  differ*  nt  azimuths  and  incidence  anyles.  Figure  ill-5  shows  records  for 
6  events  that  cover  all  azimuths.  The  S-waves  from  tin*  sedimentary  basement  vanish  for  the 
event  coming  from  an  azimuth  of  31-1  and  a  distance  of  H‘>",  indicating  normal  incidence  of  the 
P-wave  at  the  boundary.  An  89"  distant  event  emerges  at  about  18"  from  the  normal  in  a  fi-km/ 
sec  medium;  thus,  the  boundary  dips  to  the  southeast  beneath  subarray  1)2,  Other  events  with 
approximately  the  same  distance  and  azimuths  show  the  same  ext  action  of  the  S-wave,  leading 
to  the  same  boundary  configuration.  Such  a  geometry  is  in  agreement  with  the  general  trend 
determined  by  a  recent  refraction  survey 

Several  other  interesting  observations  can  be  made  fron  these  data.  Figure  lli-6  shows 
the  event  from  Az  150  .and  A  58°  with  the  SV  -component  shifted  1.5  sec  to  the  left  so  that 
the  S-wave  scattered  from  the  sedimentary  boundary  lies  directly  beneath  the  exciting  I’-wave, 

P,  l>P,  and  sP  all  clearly  show  a  corresponding  scattered  S.  The  IASA  bulletin,  based  on  just 
vertical  recordings  at  the  array,  listed  the  sP  as  pi*.  (WWSSN  data  agree  with  the  {-component 
data.,  The  pi*  phase  was  chosen  on  the  basis  of  ( I )  a  consistent  si*  time,  and  (2)  a  corresponding 
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scatterod  S-wave  for  both  pi’  and  si1  with  appropriate  amplitude  ratios  (the  S  from  pl>  should 
have  the  same  amplitude  ratio  to  pi*  as  the  S-to-I’  amplitude  ratio,  while  the  S  from  si’  should 
have  a  somewhat  larger  ratio  since  si’  is  incident  at  a  larger  angle).  l’(,  marks  the  arrival  whose 
time  and  corresponding  lack  of  a  basement  S-wave  indicate  that  it  is  probably  a  crustal  reverbera¬ 
tion  rather  than  a  part  of  the  incoming  source  waveform.  |,S^|  marks  the  S-wave  that  was  sea 
to  red  at  the  Moho,  It  conies  in  about  8  sec  after  1’.  The  same  arrival  can  he  confirmed  on  the 
various  traces  tit  Fig.  111-8  by  examining  the  S-motion  near  8  sec  after  1’  (or  I’cl’,  or  any  pri¬ 
mary  P  motion).  Sacli  a  time  delay  indicates  a  depth  between  60  and  70km  for  the  Moho  beneath 
1)2.  Similar  consistent  effects  can  he  seen  from  the  other  events  in  Fig.  111-6. 

In  summitry,  a  brief  analysis  of  5-component  seismic  motion  has  shown  promise  in  helping 
to  delineate  local  structure  and  in  separating  the  effects  of  the  structure  irom  the  incoming  source 
waveform.  As  the  itr.icture  probably  varies  considerably  over  the  array,  knowledge  of  horizon¬ 
tal  motion  at  tile  other  subarrays  would  prove  extremely  valuable  in  tm  ng  the  array  so  as  to 
achieve  its  maximum  gain.  As  only  a  single  station  experiment,  knowledge  of  the  S- waves  scat¬ 
tered  at  major  discontinuities  enhances  the  azimuth  determination  and  phase  identilirution  capa¬ 
bilities  of  the  array.  p,  lenders 

C.  LASA  SHORT-PERIOD  INSTRUMENT  RESPONSES 

To  determine  time-domain  source  characteristics  from  a  seismogram,  it  is  often  necessary 
to  remove  the  effect  of  the  seismometer  and  examine  the  motion  of  the  ground  at  the  seismometer 
site.  I  or  example,  I  easier 17  was  able  to  show  that  pP  phases  from  a  number  <>r  NTS  explosions 
were  observable  at  NOKSAU  by  removing  the  theoretical  impulse  response  from  the  data,  la  this 
SAT S,  Fils. in  and  l  rasier  have  been  ihlr  to  remove  the  responses  from  a  number  of  WWSSN  rec¬ 
ords  and  determine  the  direction  of  propagation  of  rapture  of  a  central  Asian  earthquake. 

At  1  ASA.  differences  in  signal  characteristics  are  always  observed  Irom  site  to  site,  even 
within  a  siihurray.  Some  of  these  differences  are  probable  caused  by  scattering  of  the  incoming 
seismic  signal,  while  the  remainder  are  caused  by  convolution  of  the  signal  with  different  dis¬ 
placement  impulse  responses  of  the  instruments.  To  be  able  to  separate  these  two  effects  and 
obtain  a  more  truthful  picture  of  the  incoming  signal,  it  is  obviously  necessary  to  be  able  to  re¬ 
move  the  individual  responses  from  the  data.  Since  the  response  of  eaeli  instrument  will  prob¬ 
able  change  slightly  with  tune,  an  up-to-date  measurement  of  inch  response  will  In  roipured 
whenever  the  data  need  to  he  deconvolved. 

1,1  i  <>(,(,}  we  developed  a  broadband  calibration  system  for  the  short -period  instruments 
which  has  never  been  used  to  obtain  the  response  of  a  single  instrument  to  an  impulse  in  displace¬ 
ment,  As  shown  here,  this  is  done  easily  and  can  be  done  at  any  time  to  obtain  tlm  present  im¬ 
pulse  responses  of  all  the  short-period  se, milometers, 

The  calibration  consists  of  applying  a  known  pseudo- random  sequence  of  voltage  steps  to  tin 
c  il  ibr  it  ion  coil  of  the  seismometer  and  measuring  the  simultaneous  seismometer  output  if  the 
rise  time  for  the  development  of  the  magnetic  field  produced  by  the  calibration  coil  is  short  coin- 
par  <1  with  0, nS  see  (the  sampling  rale  of  the  data),  a  step  In  voltage  applied  lo  the  coil  is  equiva¬ 
lent  to  a  step  m  acceleration  applied  m  tin*  seismometer.  In  other  words.  If  tin-  known  Input 
vol  a  .  to  the  calibration  coil  Is  \  itl,  and  the  measured  oinpiii  m  the  seismomet«  r  is  Sit)  the 
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relationship  between  S(t|  ami  V(ll  is 


S(t) 


V  (t)  l'(i  -  t)  dr 


n  ii- 1 1 


where  l,’(l)  is  the  response  of  the  seismometer  to  an  impulse  in  ground  acceleration  The  first 

derivative  of  K  -  dK/tlt  -  Is  the  seismometer  response  to  an  impulse  in  (•round  veloeltv,  while 

2  2 

Hie  seeond  derivative  of  K  -  il  I’/ilt  —  is  the  seismometer  response  to  tin  impulse  in  ground  dis- 

2  2 

placement.  Therefore,  to  obtain  d  l-'/tlt  knowing  V(t)  and  Sill,  it  is  only  necessary  to  solve 
Kip  ( 1 1 1  -  1 1  for  l’(t  I  and  to  differentiate  it  twice. 

A  calibration  of  snbarray  II)  was  made  on  2d  .fatiuary  to'2,  The  pseudo-random  voltage 
input  was  applied  simultaneously  to  all  the  It)  slnirt -period  iiutruments  for  ulmiit  S  minutes  be¬ 
ginning  at  21:07:00,  the  input  voltage  and  seismometer  outputs  being  recorded  on  tape  lig- 
ure  1 1 1  -  ~  f :» -  <  )  is  a  sample  of  the  voltage  input  and  out  put  recorded  1>^  seismometer  t  A  le,.**t- 

scpisres  approximation  to  Kit)  was  iditalned  for  each  of  the  1  f>  instruments  ami  the  subarruv  si  in 

t  o 

using  Levinson's  method  and  solving  the  eipiat ions 
M 

^  'V’vv*"  "  Tl  ' ',\sv,T>  for  v  *  0,  t. .  . . ,  M  (III-2I 

n  =  0 


where 


N-r 

-  \.._V 


t  *T  t 


<•’  (t|  *  — n - 

vv  N-r 


-,«v<Tl  ’ 


N-r 


N  -  r 


\  r  SI  1 
M  =  l‘» 


Some  of  the  Impulse  responses  I’ll)  are  shown  in  l-’ij*.  Ill-h.  while  the  result  of  rnnvol vine  the  l  it) 
obtained  for  seismometer  I  with  the  input  voltage  is  shown  In  I'ig.  Ill-7(n-r).  In  no  rase  was  the 
normall/ed  rms  difference  between  the  measured  output  and  the  result  of  convolving  the  computed 
response  with  the  input  greater  than  0.97, 

Differentiation  of  each  of  the  17  acceleration  impulse  responses  was  carried  out  in  the  fre¬ 
quency  domain  using  the  following  relations: 

f(Vl  *  \  I  ill  e*J“l  lit 

**  m 


Mil  * 

d  I  Kill  |  e 
dt 


n  L  (l' 

l  j-Cf(«,’|  ('*“’*  (to.' 
^  "  I'.f. 
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The  velocity  end  displacement  impulse  responses  «»f  some  «<f  lln'W  in'" runieiits  are  shown  in 
Kig  (tl-H  |*o  r  seismometer  I ,  the  .iitiphlod  and  phase  spectra  for  all  three  tmpiiliti  resiv.nsws 
are  shown  m  l-'tK ,  lll-,*l;i -«ll . 

I  ,0111  1  11:  Ill*h,  it  can  he  Keen  e.iHliy  that  the  res|smses  of  the  different  KelMtonim 
more  Kimll.ir  at  hlj{h  tlwui  at  low  frequencies  the  ditferences  in  the  acceleration  Impels-  . . 
spouses  are  more  obvious  titan  timse  in  displacement  tniMilse  responses.  I  he  major  differences 
appear  to  lie  in  j*n  111 :  all  these  responses  h«vc  heen  drawn  it?  he  the  Home  wire;  the  true  relative 
amplitudes  are  written  hv  the  Hltle  of  each  rcsismse  The  amplifier  at  the  «mtpnt  of  each  KeiHtnom* 
eter  in  moM  HcnHltlve  to  low -frequency  ch.in^en.  It  would  seem,  therefore,  tied  the  c.tKicHl  way 
10  check  the  adjustment  of  amplifier  ({alii  and  HelKinonieter  itimpiiiB  woultl  lie  to  ol.t  itn  I  (tl,  the 
acceleration  inipulHe  response,  hv  solving  Kq.  IIH*2|  for  each  calibration.  ^ 

Since  It  Is  now  possible  to  perform  a  Himllar  calihratlon  of  tin-  lonK-period  nHtrumentH, 

I  recommend  tl.  it  a  full  calihratlon  he  made  of  all  the  HoiHnmmeterM  at  lo\S.*  once  a  month  and 
licit  the  results  lie  available  on  tape  for  anyone  requiring  to  know  the  responses  for  that  month 

A.  Ztolkowskt 
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Fig.  lll-l(o-d).  Anolyii*  of  P-wave  oiimoth  and  slownett  onomolies  within  B-rlng  ot  LASA 
for  26  February  1967  Ewtern  Kazakh  event. 
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«-COO»DiNMC  l»»>l 

Fig.  111-2.  Graphic  disploy  of  P-wove  ozimuth  deviotions  ond  slowness 
corrections  within  C-ring  at  USA  for  26  Februory  1967  Eostern  Kazokh 


event. 
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AVERAGE  •  3  89  dB 


TRUE  AZIMUTH  (d«g) 

(c) 


Fig.  Ill-3(a-c).  Azimuth  deviation,  slowness  correction  and  relative  power  vs  true  azimuth 
for  P-wave  signals  of  31  events  measured  within  C-ring  at  LASA. 


Section  III 


Fig.  111-4.  (Vertical,  (N)arth,  and  (E)ast,  and  rotated  SV  and  SH  short-period 
seismograms  re  lorded  at  subarray  D2  showing  P  and  PcP  phases  and  corresponding 
crustal  scatter sd  S-waves  from  North  Atlantic  earthquake. 
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Fig.  MI-5.  3-component  dota  from  six  events 
from  various  azimuths  and  distances.  Note 
extinction  of  S-wove  scottered  from  sedimen¬ 
tary  section  -  tectonic  basement  boundary  for 
event  from  A z  =  314°  and  A  =  89°. 
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Section  III 


Fia  111-6  P-motion  ond  shifted  S-motion  from  event  from  Az  -  150 
ond  =  58°  showing  arrivals  of  P,  pP,  ond  sP,  and  their  correspond¬ 
ing  basement  boundary  scattered  S-woves.  Pr  ?*  crustol  reverfeerohon, 
pSM  is  P  scottered  S-wave  from  Moho. 


to) 


(b) 


[UlM 


Fiq  111-7  (o)  Pseudo-rondom  voltage  applied  to  calibration  co.l  V(t); 

(b)  ‘observed  output  ot  seismometer  1;  and  (c)  computed  output,  convo¬ 
lution  of  F(t)  with  V(t). 
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Fig.  11 1-6.  Impulse  responses  of  some  sho«t-p»'iod  seismometers  o»  S?  on  24  Jomrory  IW, 
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Fig.  III-9.  Amplitude  and  phase  spectra  for  B3 
seismometer  1  (horizontal  scale  is  frequency  in 
Hertz),  (a)  Amplitude  spectrum  of  F(t),  accel¬ 
eration  impulse  response;  (b)  amplitude  spectrum 
of  d  ( F(t)]/dt,  velocity  impulse  response;  (c)  am¬ 
plitude  spectrum  of  d^  [F(f)l/dt2,  displacement 
impulse  response;  and  (d)  phase  spectrum  of  dis¬ 
placement  impulse  response.  Phase  spectra  for 
velocity  and  acceleration  responses  have  same 
shape  but  differ  by  -ir/2  and  -ir,  respectively. 
All  vertical  and  horizontal  scales  are  linear. 


IV  EVALUATION  OF  NORSAR 


A.  SHORT-PERIOD  AMPLITUDE  VARIATIONS  AND  COHERENCY  AT  NORSAR 

The  amplitudes  of  steered  short-period  NORSAR  suharrays  have  been  measured  for  58 
events  distributed  over  16  small  regions.  Figure  1V-1  allows  the  location  of  the  regions  in 
inverse  velocity  space  relative  to  NORSAR.  The  amplitudes  measured  were  maximum  peak- 
to-peak  excursions  of  0,9-  to  3.0-llz  filte  ed  ubarray  lieams  generated  by  the  Event  Processor 
at  NOKSAK.  Only  events  with  large  signai-to-noise  ratios  (SNR)  were  used  in  order  to  avoid 
errors  introduced  by  background  noise.  The  amplitude  of  each  event  was  normalized  by  the 
geometric  mean  of  the  amplitudes  of  all  the  subarrays,  and  the  average  within  each  region  was 
calculated.  The  number  of  events  within  each  region  ranged  from  2  to  6.  The  standard  devia¬ 
tion  of  the  anomalies  within  each  region  was  also  calculated  for  each  subarray.  The  same  cal¬ 
culations  were  done  for  all  51  events  together  and  for  the  16  average  normalized  amplitudes 
from  the  different  regions.  The  results  for  four  regions  and  for  the  global  data  are  shown  in 
Table  iV-1. 

Table  IV-1  is  only  part  of  the  data,  but  it  is  representative  and  several  observations  can  be 
made.  For  a  fixed  region,  the  amplitude  variations  are  large  and  generally  quite  repea'able  for 
different  events.  However,  the  amplitude  at  a  subarray  depends  strongly  upon  the  region  —  an 
amplitude  good  for  one  region  can  be  quite  poor  for  another.  This  is  verified  by  the  global  results 
which  display  less  variation  of  average  gain  between  subarrays  and  show  larger  standard  uev 
tions  relative  to  the  average  gain. 

The  amplitude  measurements  were  made  after  new  NORSAR  station  corrections  were  intro¬ 
duced  in  early  February  1972.  It  is  worth  noting  that  the  observed  beam  amplitude,  filtered 
0.9  to  3.0  Hz,  was  on  average  equal  to  the  average  subarray  amplitude,  which  is  about  as  good 
as  could  be  expected.  An  unfiltered  beam  was  an  average  3  to  4 dll  larger  than  the  filtered 
beam,  but  the  SNR  was  worst  for  the  unfiltered  data.  Also,  one  would  expect  the  signal  loss  from 
filtering  to  be  less  for  smaller  even's  as  their  higher  frequency  content  is  more  centrally  in  the 
passband  of  the  filter. 

The  short-period  amplitude  scatter  at.  NORSAR  is  larger  than  that  at  LASA,  but  this  does 
not  mean  that  LASA  is  relatively  free  of  scatter  and  that  NORSAR  has  an  order-of-magnitude 
more  scatter.  The  true  situation  is  indicated  by  the  data  on  Figs.  IV - 2(a)  and  (b).  For  each 
subarray  at  NORSAR,  the  minimum,  maximum,  and  median  normalized  gain  overall  16  regions 
was  obtained  and  is  plotted  in  Fig.  I V - 2( a );  13  of  these  regions  are  third  zone  to  LASA.  The 
LASA  normalized  gains  for  a  single  instrument  in  each  subarray  were  available  for  these  regions.1 
Minimum,  maximum,  and  medians  were  calculated  and  are  shown  in  Fig.  IV - 2(b).  This  com¬ 
parison  of  l.ASA  and  NORSAR  scatter  is  somewhat  qualitative  and  is  between  single  instruments 
and  subarrays,  but  the  larger  scatter  at  NORSAR  is  quite  evident. 

It  is  worth  noting  that  even  subarrays  which  show  somewhat  low  amplitudes  on  average  can 
be  quite  significant  in  particular  cases.  Note  NORSAR  subarrays  4B  and  6C  in  particular.  On 
average,  they  tend  to  have  somewhat  low  amplitudes.  However,  the  maximum  gains  of  1.9  and 
2.8  shown  on  Fig.  IV-2(a)  are  for  presumed  explosions  at  the  Eastern  Kazakh  test  site.  For 
those  events,  site  6C  is  the  best  site  in  NORSAR,  and  4H  is  one  of  the  better  sites. 
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TABLE  IV-1 

AMPLITUDE  VARIATIONS  ACROSS  NORSAR 
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In  addition  to  the  amplitude  study  discussed  above,  we  have  Investigated  the  variation  in 
waveform  between  individual  NOltSAlt  subarrays  by  coni|>artng  the  coherency  of  individual  sub- 
array  beams  with  respect  to  the  full  array  beam.  Two  approaches  were  taken  within  the  coher¬ 
ency  study  a  small-population  study,  where  zero-lag  coherencies  were  computed  by  aligning 
the  subarray  beams  and  the  full  array  beam  by  eye  using  the  analysis  console;  and  a  largc- 
populatlon  study,  where  coherencies  computed  by  the  NOKSAIt  1-1  vent  Processor  were  used. 

In  the  small-population  study,  the  zero-lag  coherency  p  was  computed  using  the  short- 
period  P-wave  data  from  13  events  at  teleseismie  distances.  Once  the  i1’1  subarray  beam  x 

1.1 

and  the  full  array  beam  yj  have  been  aligned  by  ey e,p  Is  defined  as 


,  o  Vi 


Pj  /*  2  !  2\1/2 

l  1  Xii  “  yi  ) 
'10  ,J  i-0  1  / 


where  I  is  equivalent  to  5  sec  of  data.  1  he  average  at  the  j^1  subarray  was  computed  over 
the  1}  events  and  the  results  are  plotted  as  a  solid  line  in  Pig.  IY-3  where  it  is  seen  that  the 
average  value  of  p^  is  fairly  stable,  fluctuating  between  0.75  and  0.85.  Thus,  in  spite  of  the 
fact  that  certain  subarray  beams  of  certain  events  may  demonstrate  a  remarkably  low  coherence 
(p  <  0.3)  with  respect  to  the  full  array  beam,  based  on  this  data  set,  the  subarray  beams  are 
about  equally  coherent  on  the  average. 

In  the  large-population  study,  use  was  made  of  tiie  NOKSAIt  event  tape  covering  the  period 
15  February  to  15  March  1971.  Among  other  parameters  automatically  measured  from  the 
short-period  waveforms  by  the  l.'vent  Processor  are  the  coherencies  between  subarray  beams 
and  a  full  array  beam.  These  are  maximum  coherencies  on  which  the  delays  used  in  machine 
locations  are  based.  Although  they  are  not  directly  comparable  to  the  p  of  the  previous  experi¬ 
ment,  they  do  represent  a  large  body  of  data  concerning  the  coherency  of  short-per iod  wave¬ 
forms  on  the  subarray  beams.  To  elaborate,  many  of  the  events  on  the  event  tape  are  small 
with  small  SNIts;  thus,  we  should  expect  a  lower  average  coherency  at  a  given  subarray  if  the 
average  Is  based  on  a  large,  comprehensive  population.  Uc  see  that  tills  Is  the  case  in  Fig  IV- 3 
where  plotted  as  a  dashed  line  are  the  average  coherencies  of  the  individual  subarray  beams 
with  respect  to  the  full  array  beam,  the  average  being  taken  over  4 nf,  processed  "events" 
occurring  during  the  1-month  time  period.  Again  we  must  conclude,  with  one  exception,  that 
on  the  average  the  subarray  coherencies  are  quite  stable.  The  one  exception  is  so  remarkably 
low  compared  with  its  neighbors,  we  might  suspect  it  was  deleted  from  the  coherency  measure¬ 
ments  for  some  extraneous  reason  during  this  period.  It  does  not  appear  as  anomalous  in  the 
small-population  study. 

No  geographic  restriction  was  placed  on  the  406  events  used  above.  When  we  restricted 
the  events  of  the  large  population  to  0°  to  120“  azimuth  from  NOKSAIt,  the  average  coherency 
results  remained  essentially  unchanged.  Thus,  to  summarize  the  main  point,  although  certain 
subarrays  may  be  poorly  correlated  with  the  beams  for  a  given  event,  they  all  appear  to  yield 
roughly  equivalent  correlation  values  on  the  average. 


It.  T.  I.acoss 
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TABLE  IV-2 

EVENTS  L'SED  IN  NORSAR  STUDY 
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Section  IV 


0.  EVALUATION  OF  SHORT- PERIOD  NORSAR  DATA 

Hod'  wave  magnitudes  of  available  ventral  Asian  events  in  i‘>70,  recorded  by  the  interim 
Ni )HSAH  array,  have  been  sttrlied  as  jvi rt  of  the  XOHSAK  evaluation  program.  The  interim 
array  consisted  of  simile  sensors  from  most  of  the  XOHSAK  sularrays,  and  the  data  were 
diulti/cd  at  1011/.  A  simile  presumed  explosion  in  May  1’*71  is  also  included.  This  event  was 
recorded  by  the  full  XOHSAK  of  112  seismometers,  with  u  sampling  rate  of  2011/  These  full- 
array  data  are  now  hemming  available  and  will  provide  a  more  complete,  better  quality  data 
base  for  t  ie  evaluation  of  N'OUSAH 

A  list  of  the  events  is  given  in  Table  1V-2,  and  the  event  locations  are  displayed  in  the  map 
of  l  ig.  1\  -1  lour  of  the  events  are  presumed  Soviet  explosions  from  the  Soviet  test  site  near 
S«  ni|«ilatinsk  In  Eastern  Kazakh 

Array  beams  were  formed  for  each  event  bv  aligning  the  first  motions  by  eye.  The  observed 
peak-to-peak  digital  amplitude  of  each  beam  was  corrected  for  the  XOItSAU  recording  system 
and  converted  to  magnitude,  doth  the  computed  and  XOAA  magnitudes  for  each  event  are  given 
in  Table  1\  -2. 

In  Fig.  IV-5,  the  difference  between  XOAA  and  XOHSAK  magnitudes  is  shown  as  a  function 
of  epicentral  distance  from  XOHSAK  Figure  lV-fi  contains  the  same  magnitude  differences 
plotted  vs  the  period  of  I’-wave  beam.  Although  the  data  base  is  too  small  for  statistical  analysis, 
the  impression  Is  that  XOHSAK  magnitudes  are  about  0.  5  unit  lower  titan  XOAA  magnitudes, 
independent  of  the  period  and  distance  of  the  central  Asian  event.  Tills  seems  to  he  verified 
hv  event  detections  from  all  azimuths  published  In  the  weekly  XOHSAK  Seismic  liven t  Summary. 

Oil  the  other  hand.  I.acoss^  has  shown  that  I.ASA  body  wave  magnitudes  are  equal  or  slightly 
larger  than  XOAA  magnitudes,  on  the  average.  Such  XOHSAK  ami  I.ASA  magnitude  differences 
can  he  partially  accounted  for  liy  station  corrections.  A  I’-wave  arriving  at  the  surface  has  an 
energy  density  flux  which  is  proportional  to  polish),  where  p  and  o  are  the  density  and  enm- 
pressional  velocity,  respectively,  of  the  surface  rock,  and  tilt)  is  the  ground  velocity  recorded 

by  the  seismometer.  For  a  T-wavi  of  given  energy  <i.c.,  magnitude),  the  velocity  Is  proportional 

- 1  2 

to  (po)  at  the  station  site.  Thus,  hard  rock  sites,  such  as  XOHSAK,  should  measure  lower 
amplitudes  (hence  magnitudes)  than  sites  such  as  I.ASA.  which  is  situated  on  low  velocity,  sedi¬ 
mentary  rock.  1  sing  typical  values  of  p  and  o  at  XOHSAK  and  I.ASA,  we  find  that  I.ASA 
amplitudes  should  he  about  1.7  times  larger  titan  XOHSAK  amplitudes.  This  would  cause  I.ASA 
magnitudes  to  be  at  least  n.2  unit  larger  than  XOHSAK  magnitudes.  Such  an  effect  is  observed 

in  the  limited  data  shown.  . 

t  .  I' raster 

C.  VARIATIONS  OF  FALSE-ALARM  RATES  AT  NORSAR 

The  basic  detection  algorithm  is  the  same  at  I.ASA  and  NOKSAK.  The  steps  are:  (1)  I. Inear 
predetection  filtering  to  isolate  a  band  with  good  SXK;  (2)  full  wave  rectification;  (  3)  calculation 
of  a  short-term  average  (STA)  and  long-term  average  (I.TA)  by  averaging  rectified  outputs  for 
short  (typically  1.5  sec)  and  long  (typically  30  n<*v)  periods  of  time;  and,  finally,  (4)  declare  a 
detection  if  STA  l.TA  exceeds  a  preselected  threshold.  The  process  can  he  applied  at  sensor, 
subarray,  or  array  beam  levels.  It  is  clear  that  the  detection  capability  of  such  a  system  must 
decrease  if  the  noise  level  increases.  Also,  if  he  noise  level  increases  with  no  change  of 
spectrum,  the  falsc-alarm  rate  will  be  unchanged  as  long  as  the  detection  threshold  for  STA/l.TA 
is  unchanged. 
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Section  IV 


hi  fact,  the  false-alarm  rate  at  NOUS  AH  does  vary  considerably  with  background  noise 
level,  which  is  seen  quite  clearly  in  Figs.  lV-7(a)  and  (I,).  The  I.TA  is  quite  stable  and  repre¬ 
sentative  of  noise  level  at  the  output  of  the  prerectification  filter.  The  variation  in  I.TA  over 
two  time  periods  with  slightly  different  prerectification  filters  is  shown.  Intervals  contaminated 
by  obvious  events  have  been  removed.  Also  shown  is  the  number  of  detections  per  day,  the 
number  of  events  on  the  NOitSAIt  bulletin,  and  the  difference  of  these  two.  In  general,’  the  latter 
should  be  a  reasonable  estimate  <»«  false  alarms  per  .lay,  although  some  real  events  which  were 
never  studied  by  an  analyst  might  be  included.  The  detection  threshold  for  both  time  intervals 
was  the  same,  but  note  that  the  false-alar. n  rate  is  significantly  less  for  the  higher-frequency 
prerectification  filter.  Its  worst  conditions  correspond  roughly  to  the  best  conditions  of  the 
lower-frequency  filter 

An  idealized  conceptual  model  has  been  formulated  to  help  understand  these  false-alarm 
fluctuations  and  to  predict  certain  trade-offs.  First,  it  is  assumed  that  the  output  of  the  pre¬ 
rectification  filter  has  a  bandpass  spectrum  characterized  by  a  bandwidth,  center  frequei  ey,  and 
mean  square  noise  level.  It  is  also  assumed  that  the  short-term  fluctuations  of  STA/l.TA  are 
due  to  those  of  the  Sl'A,  and  that  the  I.TA  can  be  assumed  equal  to  the  average  ST  A  value.  The 
calculation  of  the  STA  from  the  rectifier  output  is  conceived  of  as  a  simple  linear  filtering  opera¬ 
tion.  It  has  been  impossible  to  obtain  the  probability  distribution  of  the  STA,  but  its  stability 
can  be  characterized  by  the  ratio  of  the  square  of  its  mean  to  its  variance.  For  a  given  threshold 
it  is  reasonable  to  expect  that  the  false-alarm  rate  will  be  a  monotone  increasing  function  of 
this  ratio.  Thus,  by  noting  the  changes  in  this  ratio  for  different  prerectification  spectra  and 

postrectification  fillers*  il  is  possible  to  predict  if  false-alarm  rates  will  increase  or  decrease. 

I  inally,  idealized  spectra  of  seismic  noise  and  gain  functions  of  prerectification  filters  have 
been  used  to  get  the  bandwidth  and  center  frequency  of  the  rectifier  input  under  various  noise 
conditions,  and  to  predict  combinations  of  seismic  noise  level  and  prerectification  filter  low- 
frequency  cutoff  which  will  result  in  equivalent  false-alarm  rates. 

At  least  for  short  periods  of  time,  the  bandpass  input  to  the  rectifier  can  be  thought  of  as 
an  amplitude-modulated  sine  wave  Alt)  sin  <2 ay  ♦  O),  where  ffl  is  the  center  frequency,  O  is 
a  random  phase,  and  Ait)  is  a  low-pass  process  with  bandwidth  H  equal  to  that  of  the  rectifier 

input.  Ait)  is  also  assumed  to  be  non  negative,  so  its  spectrum  contains  an  impulse  at  f  n 

with  area  equal  to  the  square  of  the  average  value  of  Ait).  Figure  IV-8  shows  the  power  spectral 
density  of  the  rectifier  output  under  these  assumptions.  The  size  of  the  sidebands  go  as  l/dnZ  -  l)2, 
where  n  is  the  sidebar.)  number.  The  portion  between  -ffl  and  tf  is  just  the  spectrum  of  Ait) 
multiplied  by  4/ji  .  If  H/2  >  fQ,  the  picture  is  slightly  changed  due  lo  pseudo -aliasing  since,  in 
fact,  the  rectified  spectrum  is  an  infinite  weighted  sum  of  the  spectrum  of  Ait)  shifted  by  2nf 

I  he  envelope  Ait)  would  be  reproduced  exactly,  within  a  factor  of  2/ it.  by  filtering  the 
rectifier  output  with  a  perfect  low-pass  filter  passing  frequencies  from  -f  to  4  f  The  mean 
square  value  of  this,  using  Pars  oval's  theorem,  is  (a2  .  b2|’)4/ir2  and  the  mean  value  is  a. 

Thus,  ^the^ stability  as  measured  by  the  ratio  of  the  square  of  the  mean  to  the  variance  is 

r  a  /b  I).  Hut  the  envelope  may  be  considered  to  be  Kayleigh  distributed,  for  which  we  have 
r  ir/(4  -  w).  Thus,  a2/b2 


I  7t/(4  - 


i  i. 


The  actual  filter  used  to  oh  ain  STA  from  the  rectifier  is  not  Ideal  from  -  f0  to  t  f 


It  is 


i.  >t  even  perfect  low  pass,  but  assume  for  the  moment  that  it  is  and  has  tl  e  passband  shown  in 
Fig.  IV  -H.  Taking  the  product  of  the  rectifier  output  spectrum  and  the  square  of  the  STA  filter, 
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2  ?  ,  2 

we  again  can  find  the  mean  square  value  of  ST  A  using  Parseval's  theorem  -  it  is  |(a  t  h  ).  T|l/r 

2  2  2  2  2 

The  stability  measure  is  then  rf  a  T/b  .  Substituting  for  a  'b  in  terms  of  the  bandwidth  of 
the  input  to  the  rectifier  (jives 


for  this  ideal  case  when  the  |>ostreciification  smoothing  does  not  pick  up  any  sideltand  enc/gy 
which  would  reduce  the  stat)ility.  Note  that,  all  other  factors  being  equal,  the  sideband  effect 
will  be  less  severe  for  large  fQ. 

To  complete  our  model,  it  is  necessary  to  consider  the  effective  bandwidth  and  center 
frequency  at  the  rectifier  Input  for  typical  seismic  noise  and  prerectification  filters  used  in 
N'OHSAIl.  Figure  1V-9  shows  three  power  spectral  densities  for  noise  at  NOKSAK.  They 
correspond  to  noisy,  typical,  and  quiet  conditions.  No  correction  has  been  made  for  instrument 
response.  If  !•'  is  the  log  of  frequency  and  P  is  the  number  of  decibels  down  from  the  most 
extreme  case  to  some  spectral  peak,  we  see  that  (approximately)  P  -(Ft  0.7)100.  ()n  lop. 

paper,  straight  lines  from  this  point  to  (-1.2,  -65)  and  (1.0,  -65)  give  reasonable*  approxima¬ 
tions  to  the  actual  seismic  spectrum.  The  prerectlficalion  filters  in  use  are  l>andpass  but,  for 
present  purposes,  considering  the  ranges  of  all  relevant  variables  wc  characterize  them  as 
high-pass  filters  with  rolloff  60d!l  per  decade  and  corner  frequency  f^.  Figure  IV-10  shows  the 
idealized  out|xit  spectra  from  these  filters  for  different  noise  conditions. 

The  center  frequency  and  tiandpass  of  ideal  filtered  spectra  can  be  obtained  by  simple 
geometry.  Consider  the  10-dll  bandwidth  and  center  freqjeney.  In  that  case. 


II 


(IV -2) 


fcC0,,,» 


(IV- 3) 


where 


cM(P)  =  log*1 10CIP)-  log*1 ,  rSm") 
c0(P)  Hog*1 10C(P)  ♦  log-1 1  rS-uql/z 


and 


C„»)  iH-LiI2L_ 

1  '  (100P  ♦  6500) 


<  1V-4) 


(IV -5) 


(IV -6) 


Figure  1V-11  is  a  p'ot  of  and  C0  vs  P.  They  are  roughly  the  same  shape,  and  C',  1.  32('0 

is  a  1  ;ast-squarc  .it  of  XC'0  to  using  sampler  of  P  at  0,  -  10,  -7.0,  and  -  30, 

1  he  false-alarm  rate  is  determined  by  H  and  ffl  and  is  independent  of  the  absolute  power 

level.  However,  if  f_  f  C  ( P)  is  constant  on  a  line  in  the  f  ,  P  plane,  then  so  is  H  since 

0  c  o  c  1  H 

1.32C0.  Figure  IV- 12  is  a  contour  plot  of  fn>  For  example  note  that,  if  fp  1.2  and  P  -5, 

we  should  have  the  same  false-alarm  rate  as  f  0.9,  P  -  -30,  That  is,  the  false-alarm  rate 

c 

under  best  noise  conditions  with  a  0,9-ltz  corner  should  be  the  same  as  for  worst  noise  condi¬ 
tions  and  a  1.2-tlz.  filter.  This  seems  to  be  corroborated  by  the  data  in  Fig.  IV- 7.  Also  note 
that  large  f0  implies  smaller  false-alarm  rate  since  sideband  leakage  will  be  reduced  and  [see 
Kq.(IV-t)]  bandwidth  will  be  increased. 
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Suppose  we  consider  fc  1.2.  1.6.  Then  the  effective  prerectification  bandwidths.  assuming 
V  "20‘  are  2,0J  an,J  2-6’M,z  according  to  Kqs.dV-2),  (1V-4),  and  (IV-6).  The  latter  may  be  a 
few  tenths  high  due  to  the  high-.re<|uency  cutoff  of  the  actual  filter  used,  but  that  will  be  ignored. 

,N  ,W‘  Wi‘h  '  1  5'  ,:q'  ,fV'1)  «ives  r|2  11  ’  a,lcl  '‘i2  =  t  l  7  for  the  two  filters,  respectively. 

Ibserved  values  on  a  beam  .luring  a  relatively  quiet  period  on  day  76  of  1972  were  10.7  and  14  7 
Such  loo.l  agreement  is  of  course  chance,  since  so  many  approximations  have  been  made;  but 

the  |jo jit  is  that  the  idealized  model,  which  does  not  consider  sidebands,  does  predict  the  correct 
order  for  stability. 

As  mentioned  earlier,  it  has  not  been  possible  to  determine  theoretically  the  probability 
distribution  of  ST  A  values,  experimentally  it  is  observed  to  be  quite  skew  and  clearly  not 
Gauss ia.  .  and  it  also  is  clearly  not  Hayleigh.  It  can  be  quite  well  fit  by  a  lognormal  distribution, 
lhat  is,  In(STA)  is  approximately  normal  with  nonzero  mean  p  and  standard  deviation  «r.  Values 
of  p  -1.67  and  u  0.  3  were  estimated  from  STA  values  of  a  coherent  beam  at  NOHSAll.  The 
sample  distribution  of  the  normalized  variable 

In  (STA)  -p 

y  a  ( IV  -  7) 

is  shown  in  Fig.  IV-13.  Also  shown  are  a  few  selected  points  on  the  (laussian  distribution  to 
indicate  the  good  agreement. 

If  STA  is  lognormal,  it  is  now  simple  to  estimate  expected  numbers  of  false  alarms  per 
day  on  each  beam.  Suppose  I.  is  the  decision  threshold  in  decibels.  That  is,  a  detection  is 
declared  if  twenty  times  the  base  ten  log  of  the  ratio  of  STA  to  its  long-term  average  is  greater 
than  1,.  I.et  pg  be  the  long-term  average  of  STA.  Then,  a  detection  is  declared  if 


0.1151. 


y  » 


In  -  P 


( rv-H) 


where  y  is  Gaussian  with  zero  mean  and  unit  standard  deviation.  Since  «  will  be  large,  we 
can  use  the  approximation 


i'rfy  >o)  ^  i-i/l_exp(— 2o2/»r)-|expl— 2a2/tr]  (IV-9) 

fsmg  1.  10.5,  p  and  a  given  above,  and  pg  111.6,  estimated  from  data,  we  obtain  the 

probability  of  false  alarm  for  an  observation  as  1.1  x  10"5.  It  has  been  verified  from  the 
sampled  correlation  function  of  STA  values,  which  are  calculated  every  0.5  sec  at  NOHSAli, 
hat  only  every  third  value  is  approximately  independent.  Thus,  there  are  about  0.57  x  10^ 
independent  STA  observations  on  each  beam  every  day  and  the  expected  number  of  false  alarms 

per  day,  per  beam  is  0.65.  This  is,  of  course,  for  noise  conditions  on  the  particular  day  (day 
76)  1972. 

In  fact,  on  lay  76  we  estimated  92  false  alarms  (see  Fig.  IV-7).  There  are  about 
3  0  beams  used,  and  with  0.65  false  alarm  per  beam  this  would  be  185  false  alarms  if  false 
alarms  were  independent  on  all  beams.  It  would  appear  tha*  we  must  assume  only  65  independent 
beams  to  obtain  agreement  with  observations  for  false  alarms.  A  check  of  STAs  during  detections 
on  day  76  should  indicate  how  many  beams,  on  average,  exceeded  the  threshold  for  each  false 
alarm,  but  it  has  not  been  possible  to  check  this. 
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Assistance  of  many  individuals  at  NOIISAK  and  use  of  NOKSAK  facilities  for  this  research 
is  gratefully  acknowledged.  ^  j  pacoss 

D.  COMPARISON  OF  SHORT-PERIOD  MICROSEISMIC  NOISE 
AT  LASA  AND  NORSAR 

We  have  found  that  the  P-wave  signals  at  NOKSAK  tend  to  have  much  of  their  energy  con¬ 
centrated  in  the  relatively  high-frequency  range  of  1.5  to  2.5!lz.  This  situation  is  quite  different 
from  that  at  LASA  where  the  P-wave  signal  energy  is  concentrated  primarily  in  the  0.6-  to  1.5-llz 
range.  The  relatively  high-frequency  content,  and  thus  shorter  wavelengths,  of  the  P-wave 
signal  at  NOKSAK  has  led  to  some  difficulties  due  to  loss  of  signal  coherence.  In  order  to 
achieve  the  same  short-period  detection  threshold  level  at  NOKSAK  as  at  LASA,  it  might  be 
possible  to  filter  the  NOKSAK  sig.  ils  so  that  their  dominant  energy  is  in  the  0.6-  to  1.5-llz 
range.  To  do  this,  it  is  necessary  to  know  both  about  the  signal  and  noise  levels  in  this  band 
for  both  LASA  and  NORSAK.  The  noise  problem  was  considered  first  and  the  results  of  the 
analysis  will  not  be  presented. 

The  frequency-wavenumber  spectra,  as  well  as  power  spectra  and  coherence,  were  measured 
for  the  short-period  mieroseismic  noise,  in  the  0.1-  to  5-llz  range,  at  the  Oyer  subarray  at 
NOKSAK,  as  well  as  at  LASA.  The  Oyer  subarray  consists  of  12  short-period  vertical  seis¬ 
mometers  located  "dthin  an  aperture  of  about  18  km.  The  high -resolution  method  was  used  in 
the  measurement  with  a  block  length  of  25  sec,  leading  to  a  frequency  resolution  of  0.04  Hz,  and 
with  36  blocks  which  yields  90  percent  confidence  limits  of  about  ±1.4dK.  In  addition,  the 
coherence  measurement  level  for  uncorrelated  noise  is  about  0.2,  using  the  coherence  (as  well 

4 

as  power  spectra)  estimation  procedures  described  previously. 

Some  typical  power  spectra  for  the  noise  at  LASA  and  NOKSAK  are  shown  in  Figs.  IV - 1 4(a) 
and  (b),  respectively.  The  result  for  NOKSAK  shown  in  Fig.  IV  —  14(b)  is  typical  of  six  measure¬ 
ments  which  were  made  every  other  month  during  a  1-year  period.  The  frequency -wavenumber 

5 

spectrum  for  the  noise  sample  given  in  Fig.  1V-1-  o)  was  presented  previously.  It  was  found 
that  at  0.2  and  0.4  IIz  the  noise  consisted  primari  y  of  mirface  waves  propagating  from  the  north¬ 
east  direction,  with  a  phase  velocity  of  about  3.5km/sec.  The  measurements  at  other  frequencies 
in  the  0.2-  to  1.0-llz  band  indicated  that  there  were  large  amounts  of  nonpropagating  noise,  rela¬ 
tive  to  the  amount  of  propagating  noise.  The  coherence  of  this  short-period  noise  sample  vs 
frequency  is  shown  in  Figs.  IV  —  1  5(a)  and  (b)  for  spatial  lags  of  about  1  and  3  km,  respectively. 

The  coherence  data  in  Fig.  JV-15(a)  show  that  the  1-  to  2-sec  mieroseisms  consist  of  a  nonprop¬ 
agating  as  well  as  a  propa;  ating  component.  This  propagating  component  of  the  noise  has  u 
frequency-wavenumber  spectrum  whose  structure  has  been  found  to  be  highly  diffuse  in  wave- 
number  space,  as  indicated  by  the  coherence  data  in  Fig.  IV- 1 5(b).  These  results  at  NOKSAK 
are  also  typical  of  those  found  at  LASA  for  short-period  noise. 

The  data  in  Fig.  IV- 15(a)  show  that  the  coherence,  for  a  spatial  lag  of  1  km,  increases  from 
1.5  to  4.0  Hz.  In  addition,  the  results  of  Fig.  IV - 1 5(b)  also  indicate  that  the  coherence,  for  a 
spatial  lag  of  3km,  increases  in  this  frequency  range,  but  not  to  as  large  a  value  as  that  for 
1  km.  These  data  indicate  that  there  is  propagating  noise  In  this  frequency  range  and  that  the 
relative  amount  of  such  noise  increases  as  the  frequency  increases.  In  addition,  the  data  show 
that  this  moderate  coherence  at  high  frequencies  is  definitely  due  to  ground  motion  and  not  to 
other  effects  such  as  system  noise. 
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The  power  spectra  in  Figs.  IV-14(a)  and  (b)  show  that  the  noise  levels  at  LASA  and  NORSAR, 
at  1  Hz,  are  approximately  1  and  8dB  relative  to  1  mp^/Hz  at  1  Hz,  respectively.  This  implies 
that  there  is  more  noise  energy  at  NORSAR  'it  1  Hz  than  at  LASA.  A  probable  cause  for  this  is 
the  closer  proximity  of  NORSAR  to  a  coa  uline  than  LASA.  There  is  thus  a  loss  of  about  7dB, 
relative  to  LASA,  that  would  be  incu;  x-ed  if  the  low-frequency  signal  components  at  NORSAR  were 
to  be  used  for  tiie  purposes  of  P-wave  signal  detection  even  if  the  signal  spectra  were  identical. 

However,  ‘he  signal  spectra  at  NORSAR  often  increase  with  frequency  up  to  1.5  Hz.  In 
addition,  it  is  shown  by  Frasier  u»  Sec.  IV-B  that  NORSAR  signal  magnitudes  from  Eurasia  are 
smaller  than  those  at  LASA  by  ail  average  of  approximately  10  dB.  It  is  thus  clearly  undesirable, 
both  from  signal  and  noise  considerations,  to  use  a  band  centered  o;  1Hz  for  Eurasian  events. 

J.  Capon 
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Fig.  IV-1.  Loco*ion  in  inverse  velocity,  relative  to  NORSAR,  of  16  regions  used  in  short-period 
amplitude  study. 
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NORSAR  SITE  (titered  tubarray) 


AO  B1  02  03  04  Cl  C2  C3  C4  D1  D2  D3  D4  El  E2  E3  E4  FI  F2  F3  F4 

LASA  SITE  (center  sensor) 


Fig.lV-2.  Minimum,  maximum,  and  median  subarray  relative  amplitudes  for  (a)  16  regions 
measured  at  NORSAR,  and  (b)  13  af  same  regions  measured  at  LASA. 
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23  36  «3  33  OO 

DISTANCE  (dtg) 

.  IV-5.  Difference  between  NOAA  and  NORSAR  n^,/  platted  vs  epicentral  distanci 


g.  IV-6.  Difference  between  NOAA  and  NORSAR  mj,,  plotted  vs  period  of  P-wave 
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Fig.  IV-7.  Noise  level  from  full  wove  rectifier  (LTA),  number  of  bulletin  events,  ond  estimoted  number 
of  folse  olorms  ot  NORSAR  for  (o)  17  November  to  13  December  1971  with  0.9-  to  3.  5-Hz  detection 
filter,  ond  (b)  20  Februory  to  19  Morch  1972  with  1 . 2-  to  3.  2- Hz  detection  filter. 
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Fig.lV-8.  Ideolized  power  spectrol  density  of  filtered  ond  rectified  bond-limited  noise, 
and  ideal  gain  function  of  pastrectificotion  smoothing  filter. 


[l>-t-10549'| 


1.  HIGH  NOISE  SPECTRUM 

2.  TYPICAL 

3.  LOW 

A.  P  =  -  IF  +  0.7)  100,  F  =  LOG  (») 


Fig.lV-9. 


Typical  NORSAR  ihart-period  noise  spectro  and  straight— I ine  approximations  to  them. 
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F,  LOG  (fr®qu«ncy) 


Fig.  IV-10.  Ideolized  NORSAR  noise  spectro  before  ond  ofter  filtering  with  high-pass 
filter  with  corner  frequency  fc  ond  60-dB/decode  rolloff. 


FILTERED  NOISE  PROCESS  BANDWIDTH  =  B  «  fcCg  (P  ) 
FILTERED  NOISE  CENTER  FREQUENCY  =  M  fcC  IP,) 


•50  -20  -10  0  -  30  -  20  -10  0 

UNFILTERED  PEAK  POWER  (dB  rtlotiv*  to  worit  cat*) 


Fig.  IV— 1 1 .  Scale  factors  between  bandwidth  of  rectifier  input  Cg  times  filter  corner 
fc,  ond  between  center  frequency  of  rectifier  input  C0  times  fc  os  function  of  peok 
unfiltered  noise  power  P0. 
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CONTOURS  OF  fcQUAL 
CENTER  FREQUENCY 

A  10625 
B  1  1 762 
C  1.3636 
D  1.5151 
E  16666 
F  1.8101 
G  1.9696 
H  2.1211 


Fig.  IV- 12.  Contour  map  of  rectifier  input  center  frequency-  o*  function  of  unfiltered 
peak  noise  power  and  prefilter  low-frequency  cutoff. 


Fig.  I  V-l  3.  Cumulative  distribution  of  log  of  short-term  average  of  NORSAR  beam. 
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POWER  (dB  relative  to  1  m^z/Hz  at  1  Hz) 
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14.  Typical  power  spectra  far  si 
ASA  and  (b)  NORSAR. 
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V.  SIGNALS  AND  NOISE 


A.  NONPROPAGATING  NOISE  AND  ATMOSPHERIC  LOADING 

In  the  last  SATS,  I  investigated  the  possibility  of  predicting  and  removing  some  of  the  non¬ 
propagating  noise  measured  at  LASA  by  long-period  seismometers  using  the  information  provided 
by  a  nearby  microbarograph.  This  seemed  a  reasonable  thing  to  do  as  there  is  often  a  high  coher¬ 
ence  between  the  outputs  of  a  microbarograph  and  a  seismometer  at  the  same  site  at  times  when 
the  level  of  nonpropagating  noise  is  high.  The  process  of  measuring  coherence  is  a  linear  one; 
it  would  therefore  be  expected  that  in  the  time-domain  a  transfer  ..nction  could  be  computed  which 
could  be  convolved  with  the  microbarograph  output  to  produce  that  component  of  the  seismic  noise 
which  was  coherent.  Such  a  transfer  function  was  computed  for  a  particular  noise  sample  and 
did,  indeed,  predict  all  the  coherent  noise  when  convolved  with  the  microbarograph  output.  When 
this  noise  was  subtracted  from  the  observed  noise,  the  noise  level  was  substantially  reduced.  If 
the  predictable  component  of  the  noise  had  been  largely  '  ny  effects  in  the  seismometer,  one 

would  have  expected  the  same  transfer  function  convolveu  ,  „i.  he  output  from  the  same  micro¬ 
barograph  to  have  predicted  at  least  some  of  the  noise  in  the  same  seismometer  on  another  occa¬ 
sion  when  the  coherence  between  this  microbarograph-seismometer  pair  was  again  high.  That 
this  did  not  prove  to  be  the  case  was  taken  as  evidence  tha‘  buoyancy  effects  were  negligible  and 
that  most  of  the  observed  nonpropagating  seismic  noise  was  caused  by  atmospheric  deformation 
of  the  ground. 

Since  nonpropagating  noise  is  incoherent  over  distances  greater  than  a  few  kilometers  and 
does  not  correlate  from  seismometer  to  seismometer  at  LASA,  and  since  so  much  of  it  can  be 
removed  by  merely  subtracting  that  component  which  is  correlated  with  the  microbarograph  at 
the  same  site,  it  is  clearly  the  result  of  ground  deformation  caused  by  atmospheric  pressure 
changes  in  the  near  vicinity  of  the  seismometer.  Any  far-field  atmospheric  loading  effects  must 
be  quite  insignificant.  In  describing  the  loading  process,  therefore,  we  would  expect  an  elastic 
half-space  to  be  an  adequate  model  for  the  earth,  and,  since  wind  speeds  are  two  or  three  orders 
of  magnitude  less  than  the  speeds  of  seismic  waves,  negligible  error  would  be  incurred  in  con¬ 
sidering  the  elastic  deformation  to  be  static. 

If  the  surface  of  the  half-space  is  taken  to  be  the  xy-plane  and  z  is  positive  into  the  medium, 
the  vertical  displacement  uy  at  the  point  (0,  0,  z)  to  a  pressure  load  P(x,  y)  is  given  by3 


whe  re 


1^(0,  0,  z) 


k^z  \ 

+  — — jP(x,  y)  dxdy 


r  =  (x  +  y 
2 


2  2.1/2 
+  z  )  ' 


kl  = 


1  -  v 
7rE 


k 


2 


1  t  v 
2ttE 


v  -  Poisson's  ratio 
E  =  Young's  modulus 


(V-l) 


Preceding  page  blank 
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If  the  expectation  of  is  calculated  from  Eq.(V-l)  assuming  P(x,  y)  to  be  white  and  random  in 
two  dimensions,  we  find 


0UZ21 


+ 


2 

E[P2]  rdr 


(V-2) 


This  integral  diverges  at  the  upper  limit  and  at  the  lower  limit  when  z  =  0.  In  other  words,  we 

would  expect  to  obtain  infinite  noise  from  such  a  model.  Since  this  is  not  observed,  there  is 

something  wrong  with  the  model. 

4 

llaubrich  considered  the  same  model  but  with  discrete  pressure  cells,  each  of  the  same?  area 
and  with  the  same  average  pressure  load.  The  integral  becomes  a  summation  which,  again,  does 
not  converge.  This  means,  however,  that  many  pressure  cells  are  included;  those  which  are  not 
included  always  contribute  infinitely  more  to  the  expected  noise.  The  problem  is  clear  once  we 
go  into  the  wavenumber  domain. 

If  the  applied  pressure  at  the  surface  is 

cz  =  P(o,g)  cos(2^ax  4  <fia)  cos(2n^y  +  ifi  )  (V-3) 


then  the  res 


ulting  vertical  displacement  is 


5 


u 


z 


2|T-  Z2*c  [2(1  ~v)  +2*cz| 


e-2*cz 


(V  -  4) 


where  c  =  (a2  +  p2)1/2,  and  p  is  the  shear  modulus.  It  is  clear  that  whenever  e  =  0,  u  is  infi- 

z 

nite.  This  corresponds  to  the  divergence  of  the  integral  (V-2)  at  the  upper  limit  and  is  just  a 
function  of  choosing  a  half-space  as  the  model.  If  we  had  chosen  a  spherical  earth,  this  would 
not  have  happened.  Since  we  are  sure  that  the  curvature  of  the  earth  is  irrelevant  to  this  prob¬ 
lem,  we  can  filter  out  this  D.C.  effect  by  multiplying  the  right-hand  side  of  Eq.  ( V - 4)  by  [1  -  e"yc]. 
The  divergence  at  the  lower  limit  of  the  integral  (V-2)  is  caused  by  including  infinitely  high  wave- 
numbers.  This  is  equivalent  to  modeling  the  pressure  distribution  P(x,  y)  with  an  infinite  number 
of  point  forces  which  will  always  produce  infinite  displacement  at  the  surf  a  e.  Since  point  forces 
are  not  physical,  we  can  filter  out  that  effect  by  applying  a  factor  [e’7)  c]  to  the  right-hand  side 
of  Eq.  (V-4).  With  these  additional  factors,  the  integral  (V-2)  can  be  made  to  converge. 

Examination  of  Eq.  (V-4)  provides  an  explanation  for  the  failure  of  the  transfer  function,  d< 
scribed  above,  to  remain  stable  from  one  noise  sample  to  another.  At  any  pa-ticular  wavenut  i  7  , 
the  displacement  uz  is  proportional  to  the  wavelength  l/c.  If  we  make  the  firtt-order  assumption 
that,  at  a  given  frequency,  the  wavelength  is  proportional  to  the  wind  speed,  we  should  be  able  to 
find  a  stable  transfer  function  using  the  product  of  the  atmospheric  pressure  and  wind  speed. 

Since  the  instrument  responses  of  the  mierobarograph  and  anemometer  are  entirely  different, 
these  will  have  to  be  removed  before  the  product  of  their  outputs  can  be  made. 

I  am  currently  testing  *his  hypothesis. 

J  A.  Ziolkowski 

B.  RAYLEIGH-WAVE  DETECTION  CAPABILITIES  OF  LONG-PERIOD  ALPA 
AND  LAMONT  SEISMOMETERS 


A  study  of  the  surface  wave  detection  capabilities  of  individual  long-period  seismometers 
of  the  ALPA  and  NORSAR  arrays  has  been  initiated.  The  purpose  is  to  evaluate  the  detection 
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capabilities  for  sensors  installed  at  various  sites  of  the  arrays,  and  to  compare  them  with  analog 
data  obtained  by  the  high-gain  Lamont  seismometers  installed  at  the  nearby  stations  in  Fairbanks, 
Alaska  (FBK)  and  ..t  Kongsberg,  Norway  (KON).  Preliminary  results  obtained  at  the  ALPA  and 
FBK  sites  for  31  August  1971  are  summarized  below.  Twenty-four  hour  coverage  of  digital  data 
from  several  sites  of  the  ALPA  array  and  analog  film  data  for  the  FBK  site  were  available  for 
this  day. 

Figure  V-l(a-b)  illustrates  the  relative  frequency  responses  of  the  ALPA  components  and  the 
vertical  components  of  the  system  at  station  FBK.^’7  The  relative  response  of  the  ALPA  compo¬ 
nents  peak  at  about  27  sec  and  the  relative  gains  of  the  Lamont  vertical  components  peak  at  a  gain 
of  7.7  k  at  2  5  sec  for  the  low-gain  (ZLO)  component  and  at  a  gain  of  106  k  at  32  sec  for  the  high- 
gain  (Zlll)  components.  The  ZH1  response  approximates  that  of  Pomeroy,  et  al. 

As  a  guide  to  the  probable  origin  of  events  detected  on  the  ALPA  and  FBK  records,  bulletins 
issued  by  NOS  and  SAAC  were  consulted.  Origin  times,  locations,  and  body-wave  magnitudes 
are  given  by  both  bulletins.  The  NOS  bulletin  also  reports  focal  depihs.  Distances  from,  and 
approximate  surface  wave  arrival  times  at  FBK  were  calculated  for  all  events  listed  in  the  bul¬ 
letins  assuming  a  Rayleigh-wave  velocity  of  3.5km/scc.  Table  V-l  lists  arrivals  for  12  NOS 
events,  and  Table  V-2  lists  arrivals  for  32  SAAC  events.  Some  events  appear  on  both  tables. 
Examination  of  the  ALPA  and  FBK  records  showed  indications  of  at  least  11  surface  wave  arrivals 
(hereafter  referred  to  as  ALPA/FBK  events),  some  of  which  were  clearly  discernible  on  hoth 
ZLO  and  Zlll  traces  at  FBK  as  well  as  most  channels  of  ALPA.  The  set  of  ALPA/FBK  events 
includes  some  threshold  events  which  were  not  detected  on  either  of  the  FBK  records,  but  for 
which  there  was  some  indication  of  their  presence  on  one  or  more  channels  of  the  ALPA  records; 
10  of  these  11  LP  events  could  be  associated  with  events  on  the  NOS  or  SAAC  bulletins.  Table  V-3 
summarizes  the  approximate  time  of  arrival  for  the  ALPA/FBK  events  estimated  from  the  data 
and  probable  identification  with  events  listed  in  either  the  NOS  or  SAAC  bulletins.  It  should  be 
noted  that  the  list  of  events  selected  for  discussion  in  this  report  do  not  include  a  number  of  very 
marginal  detections  which  a  more  careful  examination  of  the  records  might  confirm.  In  the  case 
of  ALPA,  increased  detection  probability  through  beamforming  can  be  expected,  and  for  FBK 
some  improvement  might  be  effected  through  digital  processing. 

As  an  illustration  of  the  relative  detection  capabilities  of  the  ALPA  and  FBK  systems, 

Figs.  V-2(a)  through  (f)  present  a  series  of  traces  of  data  recorded  by  the  two  systems  for  6  of 
the  ALPA/FBK  events.  For  each  event,  two  traces  arc  shown  for  the  LLO  and  Zlll  analog  FBK 
film  data.  Th?  traces  from  selected  vertical  channels  of  the  ALPA  sites  are  also  shown.  It 
should  be  noted  that  these  traces  comprise  short  sections  of  the  recording  for  the  entire  day  and, 
as  such,  are  not  strictly  representative  of  the  general  appearance  of  the  da  when  examined  over 
longer  time  periods.  For  example,  the  long-term  fluctuations  in  background  noise,  particularly 
for  the  Zlll/FBK  instrument,  are  not  always  fairly  represented  in  the  examples  shown. 

In  Figs.  V-2(a-f),  traces  from  the  FBK  film  records  have  the  same  rel  itive  magnification 
for  all  events.  The  amplitudes  of  the  ALPA  traces  vary  fo-  different  sites,  and  the  gains  have 
heen  adjusted  for  convenience  in  presentation. 

A  comparison  of  the  detection  capabilities  of  the  A!, PA  and  FBK  systems  on  these  data  can  be 
made  by  examination  of  the  events  shown  in  Figs.  V-2(a-f).  A  brief  Jiscussion  of  the  individual 
events  is  presented  in  the  following  paragraphs. 
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TABLE  V-3 

EVENT  IDENTIFICATION  FOR  SURFACE  WAVE  DFTECTIONS  AT  ALPA  AND  FBK 

31  AUGUST  1971 

1 

ALPA/FBK 

Events 

Arrivol 

Time 

(hours) 

Identifications  and  Locations 

1 

0043 

NOS/1  —  New  Irelond 

No  I.D. 

2 

0235 

NOS/2  —  Hindu  Kush  Region 

SAAC/4  —  Afghanistan 

3 

0621 

No  I.D. 

SAAC/7  -  North  Atlantic  Oceon 

4 

0818 

No  i.D. 

SAAC/l  3  —  North  Atlantic  Oceon 

5 

1114 

NOS/5- Jujuy  Province,  Argentina 

SAAC/l  5  —  Solto  Province,  Argentina 

6 

1323 

NOS/6  -  Near  Coast  Northern  Chile 

SAAC/l  6  —  Northern  Chile 

7 

1526 

No  I.D. 

SAAC/20  —  Kurile  Islands  Region 

8 

1856 

No  I.D. 

No  I.D. 

9 

2009 

NOS/8  -  Near  Coast  Centrol  Chile 

SAAC/26  -  Near  Coost  Central  Chile 

10 

2157 

NOS/10  -  Tongo  Islonds  Region 

SAAC/29  —  Tonga  Islands  Region 

11 

2326 

NOS/1 1  —  Guerrero,  Mexico 

SAAC/31  —  Guerrero,  Mexico 

The  extended  coda  of  a  relatively  large  magnitude  event  (mb  =  5.2)  from  the  Chile  Rise  which 
arrives  at  station  FBK  at  2  335  on  30  August  obscures  several  possible  arrivals,  from  relatively 
close  distances,  of  low  magnitude  events  located  by  SAAC  during  the  first  hour  of  31  August. 

Note  SAAC/l  from  the  Aleutian  Islands  Region. 

Event  1  [Fig.  V-2(a)]  appears  as  a  pronounced  increase  in  amplitude  at  about  OP!  hours 
[about  a  factor  of  2  for  ZLO,  and  a  factor  of  4  for  ZH1  in  the  FI1K  records  (a  factor  of  4  for  the 
ALPA  records)!.  There  is  also  an  increase  in  the  dominant  period  from  about  25  to  35  sec  at 
this  time.  Event  1  has  been  identified  from  the  NOS  bulletin  as  the  arrival  from  an  m(  4.5  event 
at  a  distance  of  81.7°.  The  SAAC  bulletin  did  not  list  this  event. 

Event  2  [Fig.  V-2(b)  1  has  been  identified  from  both  the  NOS  and  SAAC  bulletins  as  a  surface 
wave  from  the  Hindu  Kush/Afghanistan  region  of  a  magnitude  4.0  to  4.3  earthquake  and  has  an 
estimated  arrival  time  of  0232  hours.  It  is  first  visible  on  the  FBK  records  at  about  0235  and 
on  the  single  ALPA  record  at  about  023  3.  Maximum  signal  amplitude  on  either  set  of  records 
occurs  between  0237  and  0238.  Large-scale  fluctuations  in  the  background,  both  in  amplitude 
and  period  (of  the  order  of  1  minute  or  more),  make  the  ZI11  records  difficult  to  interpret  for 
small-amplitude  signals. 

Event  3  [Fig.  V-2(e)  I  has  been  identified  from  the  SAAC  bulletin  as  the  arrival  from  a  3.7 
magnitude  event  in  the  North  Atlantic  Ocean  at  a  distance  of  54°.  The  event  is  detectable  on  both 
the  FBK  and  ALPA  sensors.  No  identification  is  supplied  by  NOS. 

Event  4  [Fig,  V-2(d)  ]  is  also  identified  fiom  SAAC  bulletin  as  a  magnitude  3,7  event  from  the 
North  Atlantic  Ocean  at  a  distance  of  about  41°.  NOS  does  not  list  this  event.  As  was  the  case 
fo’.  Event  3,  it  is  clearly  detectable  on  all  sensors.  Since  both  events  have  been  estimated  to 
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have  tiie  same  magnitude  and  similar  locations,  their  records  should  be  compared  for  similar 
characteristic  features. 

Event  6  [Fig.  V -2(e)  I  has  been  identified  from  NOS  and  SAAC  bulletins  as  a  4.5  to  4.8  magni¬ 
tude  event  from  Northern  Chile. 

Event  10  [Fig.  V  —2 (f) )  has  been  identified  as  arriving  from  a  4.6  magnitude  earthquake  origi¬ 
nating  in  the  Tonga  Islands  Region.  Both  NOS  and  SAAC  lulletins  list  this  event.  The  approximate 
arrival  time  was  computed  to  be  2200  hours.  This  event  occurred  at  a  normal  depth  of  33  km. 

As  a  result  of  this  preliminary  study  of  the  AI..PA  and  FUK  systems,  it  can  be  tentatively 
concluded  that  the  surface  wave  detectic  capability  of  a  single  ALPA  sensor  is  comparable  to 
that  of  the  prejent  FBK  system.  However,  it  should  be  noted  that  the  date  chosen  for  comparison 
of  the  two  systems  was  not  ideal  in  that  neither  of  the  systems  was  necessarily  operating  at  its 
best.  Further  study  is  required  before  a  firm  statement  can  be  made. 

E.  Ashley 
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Fig.  V-2(o-f).  Long-period  vertical  motion  records  for  six  surface  wave  events  detected 
at  stations  FBK  and  ALPA;  in  each  figure,  first  two  traces  are  ZLO  and  ZHI  components 
at  FBK.  Lower  traces  are  vertical  component  records  for  one  or  more  LP  ALPA  sites. 
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Fig.  V-2(a-f).  Continued. 
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VI.  NETWORK  CAPABILITY 


A.  SHORT-PERIOD  NETWORK  CAPABILITY 

A  preliminary  investigation  into  producing  a  small-magnitudc-event  bulletin  has  been 
completed.  The  aim  of  the  investigation  was  to  produce  a  seismic  events  bulletin  of  low  mag¬ 
nitude,  >  3.5  based  on  available  data  from  NOS,  the  LASA-SAAC  bulletin,  and  the  I. ASA 
detection  log.  The  time  period  from  1  to  9  January  1970  was  chosen  as  the  investigation  interval 
since  it  would  be  so  far  removed  from  the  present  as  not  to  cause  a  problem  in  obtaining  NOS 
or  WWSSN  film  chip  data. 

In  compiling  the  event  list,  the  epicenters  reported  by  NOS  were  accepted  as  locations  that 
needed  no  verification,  However,  locations  obtained  from  NOS  data  that  were  not  published  on 
the  PDE  listings,  due  to  large  residuals  or  poor  locations,  were  examined  and  verification  of 
the  event  was  required.  Epicenters  obtained  from  the  LASA-SAAC  bulletin  that  were  not  listed 
by  NOS  were  also  considered  questionable  and  therefore  in  need  of  verification.  The  third  and 
largest  source  O'  questionable  events  was  the  LASA  detection  log.  The  process  of  event  verifica¬ 
tion  meant  that  any  questionable  event  must  be  seen  by  at  least  two  of  the  stations  in  the  selected 
list  of  WWSSN  stations  used  for  this  study  (see  Table  VI-1).  Also,  the  relocation  of  the  event 
must  fall  in  near  proximity  to  the  location  given  by  the  bulletin  or  the  detection  log. 

During  this  time  period,  NOS- PDE  listings  reported  140  events,  of  which  111  were  listed 
by  the  LASA  detection  log.  Examination  of  the  NOS  unreported  events  produced  85  verified  epi¬ 
centers,  40  of  which  were  listed  on  the  LASA  detection  log.  Those  NOS  events  not  listed  in  the 
detection  log  occurred  during  LASA  down  Tme  or  were  located  in  the  core  shadow  to  LASA,  A 
few  NOS  events  that  were  not  in  the  detection  log  were  very  small  and  reported  only  by  a  few 
near-regional  stations  or  local  networks.  The  LASA-SAAC  bulletin  reported  a  total  of  151  events, 
66  of  which  were  reported  in  the  PDE  listings. 

During  the  9-day  period,  the  LASA  detection  log  listed  1208  group  detections,  that  is,  a 
group  of  steered  beams  detecting  a  signal  within  a  10-  to  15-sec  period.  In  order  to  reduce  the 
number  of  false  alarms  and  noise  detections,  group  detections  of  less  than  5  beams  were  not 
evaluated.  This  left  a  population  of  691  detections  that  were  considered  possible  events,  or 
phases  in  need  of  verification.  The  NOS  data  verified  1  51  of  these;  the  remaining  detections 
were  considered  to  be  actual  events  if  the  signal  could  be  seen  on  two  or  more  of  the  Fielected 
WWSSN  stations.  In  this  manner,  the  detection  log  added  209  events  to  the  event  list,  The  re¬ 
maining  331  unverified  detections  were  not  necessarily  false  alarms  since  over  20  percent  of 
them  turned  out  to  be  phases  such  as  PcP,  PP,  PKKP,  etc.  Several  other  detections  were  in 
fact  real  events  and  could  be  seen  on  stations  not  in  the  selected  list  of  stations. 

A  swarm  of  events  from  Yunnan  Province,  China,  occurred  during  the  time  period  of  the 
investigation.  From  viewing  film  chips  of  station  ClfG,  Chiengmai,  Thailand,  an  additional 
34  events  were  detected  based  on  signal  character  alone,  i.e.,  separation  of  Pn,  P  ‘,  Pg,  and 
Lg  phases,  and  were  added  to  the  events  list. 

During  the  9-day  period  of  the  investigation,  a  total  of  475  events  were  verified,  an  average 
of  53  per  day.  Figure  VI-1  shows  the  relative  performance  of  each  station  used  in  the  verifica¬ 
tion  of  the  final  event  list.  R  E  Nee(!ham 

R.  M,  Sheppard 


Preceding  page  blank 


83 


Section  VI 


TABLE  VI-1 

STATIONS  USED  FOR  VERIFICATION  PURPOSES 

IN  THE  NETWORK  CAPABILITY  STUDY 

WWSSN  Stations 

— 

Gain 

(in  thousands) 

Code 

Station 

Region 

SP 

LP 

AFI 

Afiamalu 

Samoa 

12.5 

0.750 

ARE 

Arequipa 

Peru 

50.0 

1.500 

ATU 

Athens  University 

Greece 

12.5 

1.500 

BUL 

Bulawayo 

So.  Rhodesia 

100.0 

1.500 

CHC 

Chiengmai 

Thailand 

200.0 

3.000 

COL 

College  Outpost 

Alaska 

100.0 

1.500 

EIL 

Eilat 

Israel 

200.0 

3.000 

KBL 

Kabul 

Afghanistan 

400.0 

6.000 

MAT 

Matsushiro 

Honshu 

100.0 

3.000 

PMG 

Port  Moresby 

New  Guinea 

50.0 

3.000 

PRE 

Pretoria 

South  Africa 

50.0 

1.500 

QUb 

Quetta 

Pakistan 

200.0 

6.000 

SJG 

San  Juan 

Puerto  Rico 

50.0 

0.750 

SPA 

South  Pole 

Antarctica 

100.0 

0.375 

TAU 

University  of  Tasmania 

Tasmania 

25.0 

0.750 

TRN 

Trinidad 

West  Indies 

25.0 

1.500 

TUC 

Tucson 

Arizona 

200.0 

1.500 

Stations  Other  Than  WWSSN  Stations 

LAO 

Lasa  Array 

Montana 

I  MBC 

Mould  Bay 

NW  Territory 

100.0 

5. 000 
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Section  VI 


B.  INTERNATIONAL  SEISMIC  MONTH 

An  effort  in  under  way  to  produce  a  very  comprehensive  seismic  bulletin  for  the  period  of 
20  February  through  19  March  1972  for  the  purpose  of  evaluating  some  of  the  large-scale  seismic 
systems  that  have  recently  become  operational.  The  cooperation  of  several  other  agencies  (see 
below)  was  enlisted  in  order  to  support  tlie  International  Seismic  Month  (ISM)  effort.  Lincoln 
Laboratory  is  acting  as  host  for  the  experiment.  The  first  phase  of  the  ISM  effort  is  to  produce 
as  complete  a  list  of  epicenters  as  is  possible  using  the  detection  logs  from  both  I.ASA  and 
NOllSAlt,  and  arrays  at  llagfors,  Gauribidanur,  Warramunga  and  Yellowknife.  'I  hese  detection 
logs,  coupled  with  analyst  reports  from  the  single  station  in  the  Canadian  network  and  selected 
WWSSN  sites,  are  providing  the  basic  data  for  the  ISM  bulletin.  Progress  on  the  bulletin  has 
been  promising  thus  far,  and  full  use  has  not  yet  been  made  of  I.ASA  and  NORSAU  detection  logs. 
Even  so,  the  number  of  events  in  the  bulletin  is  well  over  a  thousand.  We  expect  that,  upon 
completion  of  the  bulletin,  the  total  number  of  events  will  approach  two  thousand.  At  various 
stages  during  compilation  of  the  bulletin,  listings  are  being  supplied  to  those  cooperating  in  the 
study. 

The  preparation  of  the  bulletin  is  only  one  phase  of  the  ISM.  Insofar  as  possible,  detailed 
measurements  of  long-  and  short-period  data  and  potential  discriminants  will  be  accumulated  to 
form  the  complete  data  base.  All  participants  will  then  have  full  access  to  these  data  to  evaluate 
detection  and  discrimination  capabilities  of  individual  sites,  subsystems,  and  the  overall  system, 
which  can  be  conceived  of  as  the  current  deployed  capability.  Presumably,  some  of  the  sites 
used  (or  similar  ones)  could  be  incorporated  into  a  system,  if  that  were  required.  The  study 
should  also  indicate  how  current  capability  is  deficient. 

The  following  organizal  ions  are  participating: 

Department  of  Energy,  Mines,  and  Resources,  Canada 
l.amont- Doherty  Geological  Observatory 

National  Earthquake  Information  Center,  Boulder,  Colorado 
NORSAR,  Norway 

Research  Institute  of  National  Defense,  Sweden 
Texas  Instruments,  Alexandria,  Virginia 
United  Kingdom  Atomic  Energy  Authority,  England 
SAAC,  Alexandria,  Virginia 
M.l.T.  Lincoln  Laboratory 

R.  M.  Sheppard 
R.  E.  Needham 
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